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We study the generation of a non-equilibrium plasma in scalar QED with iV-charged 
scalar fields in the cases: a) of a supercooled second order phase transition through spin- 
odal instabilities and b) of parametric amplification when the order parameter oscillates 
with large amplitude around the minimum of the potential. The focus is to study the non- 
equilibrium electromagnetic properties of the plasma, such as photon production, electric 
and magnetic screening and conductivity. A novel kinetic equation is introduced to compute 
photon production far away from equilibrium in the large N limit and lowest order in the elec- 
tromagnetic coupling. During the early stages of the dynamics the photon density grows ex- 
ponentially and asymptotically the frequency distribution becomes N p h(u) ~ a m 2 /[A 2 lo 3 ] 
with A the scalar self-coupling and m the scalar mass. In the case of a phase transition, elec- 
tric and magnetic fields are correlated on distances £ (t) ~ ^Jtjm during the early stages of 
the evolution and the power spectrum is peaked at low momentum. This aspect is relevant 
for the generation of primordial magnetic fields in the Early Universe and for photopro- 
duction as a potential experimental signature of the chiral phase transition. Magnetic and 
Debye screening masses are defined out of equilibrium as generalizations of the equilibrium 
case. While the magnetic mass vanishes out of equilibrium in this abelian model, we intro- 
duce an effective time and wave-number dependent magnetic mass that reveals the different 
processes that contribute to screening and their time scales. The Debye mass turns out to 
be m 2 Dl , b ~ a m 2 /X for a supercooled phase transition while in the case of an oscillating 
order parameter an interpolating time dependent Debye mass grows as avrot/A due to a 
non-linear resonance at low momentum in the charged particle distribution. It is shown how 
the transverse electric conductivity builds up during the formation of the non-equilibrium 
plasma. Its long wavelength limit reaches a value <7fc~o ~ ti m/A at the end of the stage 
of linear instabilities. It is shown that the electric conductivity stays finite for all k in- 
cluding k = for finite time. In the asymptotic regime it attains a form analogous to the 
equilibrium case but in terms of the non-equilibrium particle distribution functions. 
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I. INTRODUCTION AND MOTIVATION 



The study of the dynamics of phenomena strongly out of equilibrium is very relevant in cosmology where 
it plays a fundamental role in the consistent description of inflationary scenarios, baryogenesis and of gener- 
ation of primordial magnetic fields. Also in Relativistic Heavy Ion collisions where it now acquires further 
phenomenological importance since the Relativistic Heavy Ion Collider (RHIC) at Brookhaven begins oper- 
ation. RHIC and the forthcoming Large Hadron Collider at CERN will probe the quark-gluon plasma and 
the chiral phase transitions in an extreme environment of high temperature and density. These experimen- 
tal programs have inspired intense theoretical efforts to understand the formation, evolution and potential 
experimental signatures of the quark-gluon plasma as well as relaxation and transport phenomena on 
unprecedented short time scales. There are several fundamental questions which define to a large extent 
the theoretical aspects of this program: how does the quark-gluon plasma form and equilibrates from the 
evolution of the parton distribution functions? what are the time scales for electric and magnetic screening 
that dress the gluons and cut-off small angle scattering? how does a hydrodynamic picture of the space-time 
evolution of the plasma emerge? what are the experimental signatures?. These and other fundamental but 
extremely difficult questions are being addressed from many different perspectives. An important approach 
that seeks to describe the space-time evolution of partons is based on transport equations that describe 
partonic cascades starting from a microscopic description and incorporate semi-phenomenologically some 
screening corrections in the scattering cross sections [f|-p| • A correct description of electric and magnetic 
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screening is very important in this program since both act as infrared cutoffs in transport cross sections and 
determine energy losses in the plasma. Amongst the several potential experimental signatures proposed to 
detect the QGP, photons and dileptons are deemed to be clean probes of the quark-gluon plasma because 
they only interact electromagnetically an d their mean-free paths are much larger than the size of the 

fireball w 20fm. Hence, these electromagnetic probes could provide clean signatures of equilibration or out 
of equilibrium phenomena unhindered by the strong interactions. Non-equilibrium phenomena associated 
with a quenched chiral phase transition could have potentially important electromagnetic signatures in the 
photon spectrum if there are strong charged pion fluctuations during the phase transition. A preliminary 
study in this direction was pursued in 0] where it was indicated that departures from equilibrium in the 
photon distribution at low momentum could provide a signature of a supercooled chiral phase transition 
]S| . In cosmology, post-inflationary phase transitions or the fast evolution of an inflaton field after inflation 
could generate the hot plasma that describes the standard big bang scenario with a radiation dominated 
Friedman- Robertson- Walker cosmology at the end of inflation Furthermore, non-equilibrium effects dur- 
ing cosmological phase transitions had been conjectured to generate the primordial magnetic fields that could 
act as seeds to be amplified by dynamo mechanisms as an explanation for the observed galactic magnetic 
fields 10 Theoretical models for generation of primordial magnetic fields involve strong fluctuations of 



charged fields that lead to non-equilibrium electromagnetic currents |T^-jlJ], much like the strong fluctua- 
tions in the pion fields during a possible supercooled chiral phase transition and the possibility of photon 
production associated with these fluctuations . 

Thus we see that physically relevant non-equilibrium physical phenomena are common to cosmology and 
the quark-gluon plasma and chiral phase transition and it has been conjectured that indeed primordial 
electromagnetic fields can be generated from strong electromagnetic fluctuations at the quark-hadron phase 
transition fll5| , [T6| . An important ingredient both in the quark-gluon plasma as well as in the formation of as- 
trophysical and cosmological plasmas is a description of the transport properties, in particular the screening 
masses and the electrical conductivity. Screening masses are an important ingredient in charmonium sup- 
pression which is one of the potential probes of the QGP jlTj and regulate the infrared behavior of transport 
coefficients 

The electrical conductivity plays an important role in the formation and correlations of primordial mag- 
netic fields in the early universe and contributes to ohmic heating and therefore energy losses and entropy 
production in the QGP. The electrical conductivity in the early universe was estimated in and (equi- 
librium) screening corrections were included in |2C| ] . More recently the electrical conductivity of the plasma 
at temperatures near the electroweak scale was calculated in ^ij including Debye and dynamical (Landau 
damping) screening of electric and magnetic interactions. 

Hence, there are common relevant problems in cosmology, astrophysics and ultrarelativistic heavy ion 
collisions that seek a deeper understanding of the physics of the formation of a plasma beginning from a 
non-equilibrium initial state of large energy density, its evolution, the onset of electric and magnetic screening 
phenomena and the generation of seeds of bulk electric and magnetic fields, i.e. photon production. 

A first principles description of the formation of a hot plasma and its dynamical evolution from an initial 
state of large energy density beginning from QCD or the Standard Model would be a desirable goal, but 
clearly an extremely complicated task. 

The goals of this work: 

In this article we study a model that bears many of the important aspects of QCD and the Standard 
Model which combined with a non-perturbative framework allows us to provide quantitative and qualitative 
answers to many of the questions associated with the formation and evolution of a non-equilibrium plasma. 

The model that we propose to study is scalar QED with A-charged scalar fields coupled to one U(l) 
photon field and one neutral scalar field that plays the role of an order parameter for a phase transition. The 
model is such that the U(l) local gauge symmetry associated with the photon field is not spontaneously 
broken much in the same manner as the usual electromagnetic field in the Standard Model. Besides, this 
model being a suitable framework to study the questions posed above, we will argue that it is potentially 
relevant to the description of photon production during the chiral phase transition of QCD. Therefore, the 
dynamics and mechanisms revealed in this model could prove to be very valuable in the description of the 
generation of primordial magnetic fields during one of the QCD phase transitions in the early universe and 
also in photon production during the chiral phase transitions in heavy ion collisions. 
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Furthermore scalar QED has been shown to share many properties of spinor QED and QCD in leading 
order in the hard thermal loop approximation |p7| , j32f , hence the model studied in this article can serve as a 
useful and relevant testing ground to study similar questions in QED and QCD. 

Since the non-equilibrium processes that lead to the formation of the plasma are non-perturbative, we 
resort to the large N limit as a consistent framework to study the non-perturbative dynamics. We take 
the electromagnetic coupling to be perturbative and compute various quantities, such as the rate of photon 
production, magnetic and Debye masses and the transverse conductivity to leading order in the large N limit 
and to lowest order in the electromagnetic coupling, discussing the validity of weak coupling in each case. 

The focus of this work centers on the following aspects: i) The description of the formation of a non- 
equilibrium plasma of charged particles during a stage of strong non-equilibrium evolution beginning from 
an initial state of large energy density, ii) The production of photons and therefore of electric and magnetic 
fields from the strong fluctuations of the charged fields. This aspect is relevant for the formation of primordial 
magnetic fields in the early universe and also for photon production during non-equilibrium stages for example 
of the chiral phase transition, where the charged fields would be the pions. iii) The dynamical aspects of 
electric and magnetic screening. We study in detail the magnetic and Debyc masses and the time scale of the 
different processes that contribute to screening, iv) The non- equilibrium transverse electrical conductivity. 
We analyze in detail the build-up of conductivity as the plasma is forming and its asymptotic limit, comparing 
to the equilibrium case. 

In particular two important situations are studied: a): a 'quenched' (or supercooled) second order phase 
transition in which the initial state of large energy density is the false vacuum (the quantum state is localized 
at the top of the potential). The dynamics in this case is described by the process of spinodal decomposition 
and phase separation, characterized by the exponential growth of long-wavelength unstable fluctuations. 
These instabilities and the ensuing large fluctuations of the charged fields and particle production result in 
the formation of a non-equilibrium plasma and the non-perturbative production of photons and therefore of 
electric and magnetic fields. The spinodal instabilities are shut-off by the non-linearities and the resulting 
plasma possesses a non-equilibrium distribution function of charged scalars peaked at low momenta, b): 
The stage of large amplitude oscillations of the order parameter around the minimum of the potential. This 
stage arises for example after a phase transition in which the order parameter has rolled down the potential 
hill and is oscillating around one of the minima of the potential. Such would be the case in the case of 
the chiral phase transition where a small explicit symmetry breaking term (that gives mass to the pions) 
will force the isoscalar order parameter to evolve towards the minimum. This stage is characterized by 
parametric amplification of quantum fluctuations of the charged fields and again results in non-perturbative 
production of charged scalars and of photons . This stage is also relevant in cosmology and describes 
the reheating process after an inflationary phase transition or in chaotic inflationary models Q. The 
phenomenon of parametric amplification of quantum fluctuations during the oscillatory phase of the order 
parameter, the inflaton in the cosmological setting, has been recognized as a very efficient mechanism of 
particle production and reheating in the early universe [p4|]25| l . Parametric amplification of pion fluctuations 
after a supercooled chiral phase transition has also been recognized to be an important possibility in heavy 
ion collisions [p6| . Both non-equilibrium phenomena are non-perturbative in the scalar quartic self-coupling. 
Therefore, the dynamics in the scalar sector is studied consistently in leading order in the large N expansion, 
while electromagnetic phenomena are studied to lowest order in a. 

Spinodal instabilities or parametric amplification of quantum fluctuations of charged fields result in the 
formation of a non-equilibrium plasma. In both cases strong fluctuations in the electromagnetic currents 
result in the production of photons i.e. electric and magnetic fields as well as screening currents generating 
screening masses and an electrical conductivity in the medium. 

Thus, our main objectives are to study the dynamics of formation of the non-equilibrium plasma, photon 
production and the power spectrum in the generated electric and magnetic fields, the onset of electric and 
magnetic screening phenomena described in real time and the build up of conductivity in the medium. 
Equilibrium aspects of hot scalar QED had been previously studied p7 
equilibrium aspects to the equilibrium case to highlight the differences and similarities. 

Results: a) Photon production: 

We have derived a consistent kinetic equation to describe photon production in situations strongly out of 
equilibrium and used this equation to lowest order in a (the electromagnetic coupling) and leading order 



28 1 and we will compare the non- 
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in the large N limit for the charged fields, to obtain the spectrum of photons produced via spinodal and 
parametric instabilities. In the case of spinodal instabilities which correspond to the case of a supercooled 
(second order) phase transition we have obtained the power spectrum and correlation function of the electric 
and magnetic fields generated during the non-equilibrium stage. We find that there is a dynamical correlation 
length that grows as £(t) ~ yi at short times. It determines the spatial correlations of the electromagnetic 
fields. The power spectrum is peaked at long- wavelength with an amplitude ~ a/ A 2 with A the quartic self- 
coupling of the charged scalar fields. In the case of parametric amplification the power spectrum peaks near 
the center of parametric resonance bands; the amplitude being also ~ a/A 2 but the electric and magnetic 
fields have small correlation lengths. In the asymptotic regime the distribution of produced photons as 
function of frequency to behaves as ~ a m 2 /[uj 3 A 2 ]. This entails a logarithmically infrared divergent number 
of photons but a finite total energy. In the case when the plasma is generated by spinodal instabilities, 
the asymptotic photon distribution continues to grow proportional to log mr due to collinear singularities. 
These behaviors points to the necessity of a resummation perhaps via the dynamical renormalization group 
introduced in reference p6| ], 

b) Magnetic and Debye screening masses: We introduce a definition of the magnetic and Debye 
screening masses out of equilibrium which are the natural extension of that in equilibrium p9| , |31 32 1 . We find 
that the magnetic mass out of equilibrium vanishes at order a through cancellations akin to those that take 
place in equilibrium. Furthermore, we introduce an effective magnetic mass that describes non-equilibrium 
screening phenomena for long-wavelength fluctuations as a function of time and which reveals the different 
time scales of the processes that contribute to the cancellation of the magnetic mass. Asymptotically for long 
times and in the long-wavelength limit we find that processes which are the non-equilibrium counterpart of 
Landau damping contribute on time scales which are much longer than typical production and annihilation 
processes. 

The extrapolation of this time dependent effective magnetic mass to the zero momentum limit at finite time 
reveals an unexpected instability in the time evolution of transverse electromagnetic mean fields during the 
time scales studied in this article. This is a rather weak instability presumably related to photon production 
although the precise relation is not clear and deserves further study. 

In the case of spinodal instabilities we find that the (electric screening) Debye mass at leading order in 
e 2 and 1/N is finite and given by m 2 Deb = 8 |mij| 2 e 2 /A + 0(1). In the case of parametric amplification 
the Debye mass grows monotonically with time as y/t times a coefficient of order 0(e 2 /A). This result is a 
consequence of non-linear resonances |37 which make the charged particle distribution strongly peaked at 
small momentum, Afk(t — oo) ~ |mfl| 2 /[A fc 2 ]. Since the Debye mass is determined by the derivative of the 
distribution function, the singularity at small momentum results in a divergent Debye mass for asymptotically 
long time. This result, valid to first order in a, strongly suggests that a resummation of electromagnetic 
corrections will be required in the case of parametric resonance. Such a program lies outside the scope of 
this work and will be the subject of a forthcoming study. 

c) Transverse electric conductivity: As the plasma of charged particles forms the medium becomes 
conducting. We study the transverse electrical conductivity from linear response out of equilibrium (Kubo's 
conductivity) as a function of time to lowest order in the electromagnetic coupling. The early time behavior 
during the stage of spinodal or parametric instabilities results in a rapid build up of the conductivity which 
attains a non-perturbative value 0(a m/A) at the end of this stage. We find that the conductivity is finite 
for all k (including k = 0) at finite time. Asymptotically at long times, the conductivity attains a form sim- 
ilar to the equilibrium case (to lowest order in a) but in terms of the non-equilibrium distribution functions. 



This feature of the asymptotic conductivity must apply to other physical magnitudes for asymptotic times. 
Namely, one can compute their t — > oo limit just replacing the thermal occupation numbers in their equilib- 
rium expression by the out-of-equilibrium distribution functions. 

The article is organized as follows: in section II the model is introduced and the large N limit is described. 
In section III we review the main features of spinodal decomposition and parametric amplification and 
introduce the relevant non-equilibrium Green's functions necessary for the calculations. In section IV we 
study photon production both during the early stages of the instabilities as well as at asymptotically long 
times, In section V we study photon production in equilibrium to contrast and compare to the non-equilibrium 
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results. In section VI we study magnetic screening and the magnetic mass out of equilibrium. Just as in the 
equilibrium case in this abelian theory, we show that the magnetic mass vanishes, but point out the different 
time scales for the processes involved. A suitably defined effective magnetic mass describes non-equilibrium 
aspects of magnetic screening on intermediate time scales. Section VII studies the Debye (electric) screening 
mass, and it is argued that in the case of parametric amplification the Debye mass diverges because as a 
result of a singular distribution function for the charged scalars at low momentum. In section VIII we study 
Kubo's (linear response) transverse electrical conductivity to lowest order in a. In particular we focus on the 
build-up of conductivity during the early stages of formation of the plasma. We compare the conductivity 
in the asymptotic time regime to the result in equilibrium. Our conclusions are summarized in section IX. 
Here we also discuss the limit of validity of our studies and the potential phcnomenological implications 
of the results from this model. An Appendix is devoted to a novel kinetic equation that describes photon 
production away from equilibrium. 



II. THE MODEL: SQED WITH N CHARGED SCALARS IN THE LARGE N LIMIT 

We focus on the non-equilibrium dynamics of the formation of relativistic quantum plasma at high density 
after a phase transition, cither via long-wavelength spinodal instabilities in the early stages of a rapid 
(quenched) second order phase transition or by parametric amplification of quantum fluctuations as the 
order parameter oscillates around the equilibrium minimum. Previous work P5| , ^6||37| revealed that both 
types of phenomena are non-perturbative in the scalar self-coupling, hence we propose to use the large N 
limit as a consistent tool to study non-equilibrium phenomena non-perturbatively. Our main goals are to 
provide a quantitative understanding of several important processes that are of interest both in cosmology 
as well as in the formation of a quark-gluon plasma: i) non-equilibrium production of photons, i.e. the 
non-equilibrium generation of electromagnetic fields, ii) the dynamics of screening and generation of electric 
and magnetic masses strongly out of equilibrium, iii) the build-up of conductivity in the non-equilibrium 
plasma. 

We consider a version of scalar quantum electrodynamics with N charged scalar fields $ r to be collectively 
referred to as pions coupled to a neutral field a is such a way that the scalar sector of the theory has 
an 0(2N + 1) isospin symmetry. The coupling to the electromagnetic field reduce this symmetry to an 
SU{N) g i bai x U(l)i oca i. When we consider the breaking of the isospin symmetry, the neutral scalar field 
a will acquire an expectation value, but not the charged fields $ r . There are two main reasons for this 
choice a) this allows to separate the Higgs phenomenon and generation of mass for the vector field from 
truly non-equilibrium effects and b) we seek to describe a phenomenologically relevant model, in particular 
the role of non-equilibrium pion fluctuations during the chiral phase transition wherein electromagnetism is 
not spontaneously broken by chiral symmetry breaking. 

The same methods can be used to study the Higgs phenomenon out of equilibrium and we expect to report 
on such study in the near future. Furthermore, as we seek to describe some relevant phenomenology for low 
energy QCD, this model describes the large N limit of the 0(4) gauged linear sigma model that describes 
the three pions. Electromagnetism is unbroken but isospin is broken by the coupling of the charged pions to 
electromagnetism and this is captured by the model under consideration. 

In this abelian theory it is straightforward to provide a gauge invariant description by requiring that the 
set of first class constraints, Ho = ; V • E — p = annihilate the physical states [p8| , p9| with Ho being 
the canonical momentum conjugate to the temporal component of the vector field and V ■ E — p = is 
Gauss' law and p is the charge density. This procedure is described in detail in |38|,|39| where it is shown 
to be equivalent to a gauge-fixed formulation in Coulomb's gauge. The instantaneous Coulomb interaction 
is traded by a Lagrange multiplier Aq (x, t) not to be confused with the original temporal component of the 
gauge field. The issue of gauge invariance is an important one because we will study the distribution function 
of charged scalar fields and by providing a gauge invariant description from the beginning we avoid potential 
ambiguities. 

In this formulation we introduce the physical fields 



6 



{a,$ r ,<S>l,X T ,Ao), r=l,...,N 



The electromagnetic potential is a physical field which satisfies the transversality condition 

V ■ A T = 0, 

whereas ylo is the Lagrange multiplier associated with the Gauss' law constraint 

V • E = -V 2 A = p. 

Thus Aq is a non-propagating field completely specified by the charge density evolution. 
To simplify expressions, we now use the following notations, 



N N N N 



r—1 r—1 



With these notations the Lagrangian density is written 

C = d + C 2 + £3 (2.1) 



with 



(2.2) 



C 2 = ^d fl A T -d"A T + ^(VA f (2.3) 



£ 3 = — ■ ($tv$ - V$ f $) (A T - A 2 ) _ ^ Q /$$t _ $t$ (2.4) 



and 



We have rescaled the couplings with the proper powers of N so that e , A are fixed in the large N limit. This 
rescaling allows a consistent identification of terms as powers of 1/N. 

We focus on the evolution of initial states with a nonperturbatively large energy density ( of order m i /X 3> 
to 4 ) in two different situations: i) m 2 < 0: this case corresponds to a symmetry breaking potential. We 
will choose the neutral field a to describe the direction of global symmetry breaking, hence the local gauge 
symmetry describing electromagnetism is not spontaneously broken, i.e. ($) = 0. A rapid (quenched 
or supercooled) symmetry breaking phase transition can be described by assuming that m 2 changes sign 
suddenly from positive describing a symmetric potential to negative describing a symmetry breaking potential 
p6| , p7|j34| . The long- wavelength modes become unstable and grow exponentially, this is the process of 
spinodal decomposition and the hallmark of phase separation. 

We emphasize that in the case under consideration the choice of negative sign does not result in the 
spontaneous breakdown of the gauge symmetry, since in this model the gauge field does not couple to the 
order parameter a. Indeed, the global gauge symmetry is always unbroken as the charged fields do not 
acquire an expectation value. 

ii) m 2 > , {cr)(t = 0) w m/vA: in this case the expectation value of the sigma field will oscillate inducing 
large parametric amplification of the <f> field. In both cases the quantum fluctuations of the fields will become 
non-perturbatively large in the scalar self-coupling and these will be treated in the leading order in the large 



N limit (mean field) [p6 37 1. The electromagnetic interaction instead, being of order a will only give small 



corrections to the scalar field evolution: thus the backreaction of the gauge field on the evolution of the 
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scalar field will be neglected. Therefore to leading order in N the equations of motion for the scalar sector 
are the same as those obtained in refs. |37| in absence of electromagnetic coupling. 

Assuming in general that the sigma field acquires a non-equilibrium expectation value we shift 

a(x,t) = VN^(t)+ X (x,t) ; ( X (x,t))=0 (2.5) 

where the expectation value is taken in the time evolved density matrix or initial state. The large N limit 
in the scalar sector can be obtained either by introducing an auxiliary field [|40|| or equivalently in leading 
order by the Hartree-like factorizations ]36| ] 

($ f $) 2 -> 2($ t $)$ t $ (2.6) 
x $t$ x ($t$) (2J) 

The non-linear terms of the a field lead to subleading contributions in the large N limit, and to leading 
order the dynamics is completely determined by the N complex scalars $. The factorization that leads 
to the leading contribution in the large N limit makes the Lagrangian for these fields quadratic (in the 
absence of the gauge coupling) at the expense of a self-consistent condition: thus charged fields <£> acquire a 
self-consistent time dependent mass. The dynamics is determined by the equation of motion of ip(t) and by 
the Heisenberg equations of the charged fields. 



III. SPINODAL AND PARAMETRIC INSTABILITIES: SUMMARY OF MAIN FEATURES 



Before we begin our study of non-equilibrium photon production and the emergence of dynamical masses, 
we review the main features associated with the non-equilibrium dynamics of the scalar fields to provide the 
physical picture and the basic ideas upon which we will elaborate with the inclusion of the gauge fields. For 
more details the reader is referred to p6p7| ] . As mentioned above the leading order in the large N limit can 
be obtained by a Hartree-like factorization that turn s the Lagrangian into a quadratic form. The equation 
of motion for the expectation value <p(t) [see eq.(2.5)] is given by 



<p{t) + m 2 tp{t) + ~ <p 3 (t) + ^($t$) tp(t) = 



(3.1) 



Introducing the usual decomposition 



d 3 k 



a r (k)h{t) e ik -* + bt(k)r k (t) e~ ik - s 



(3.2) 



*J(t,f) 



d 3 k 



b r (k)f k (t) e ik - s + at(k) f* k (t) e"* 



7 k ■ X 



(3.3) 



we find that the charged fields obey the Heisenberg equations if the mode functions /& (t) obey the following 
equations of motion p6p^] 



^ + fc 2 + m 2 + ^ 2 (<) + A <(i >t $) 



fk(t) = 



(3.4) 



We will choose the initial state to be the state annihilated by the a r (k), b r (k) operators and determined by 
the following initial conditions on the mode functions, 



A(o) = 



l 



A(0) = -iW k /fc(0) 



(3.5) 
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The frequencies Wk will be chosen in the particular cases to be analyzed below. This choice of initial state 
with the initial conditions given by (^^) corresponds to the vacuum of the Fock quanta of oscillators of 
frequencies Wk ■ This initial state can be generalized straightforwardly to a thermal density matrix, but the 
main physical mechanisms can be highlighted in a simpler manner by the choice of this state. With this 
choice one finds 



N [ ' A J (2tt 



d 3 k 
f 



(3.6) 



This expectation value is ultraviolet divergent, therefore the renormalization must be carried out consistently 
in terms of mass and coupling constant and the reader is referred to |36|j3^ 1 for details. 

It proves convenient to introduce dimensionless variables in terms of the renormalizcd mass and coupling 



\m R \ t 



\m R \ 



8tt 2 



n q = 



w k 

\m R \ 



(3.7) 



Xl! <p 2 (t) , <p q (T) = \m R \?f k (t) 



2|m fl |2 

and the subtracted self-consistent self-energy fl36| , |37| 



gZ(r) =g I q 2 dq\ K(r)| 2 - |^(0)| 2 + ^ ~ ' ' 



2q 3 



-77 2 (0)+77(t) 2 + 5 S(t) 



(3.8) 



(3.9) 



From now on we set the only dimensional variable in the problem \m R \ = 1 and all dimensionful quantities 
will be in units of |m/j|. 

To leading order in the large N limit the dynamics is completely determined by the following equations 
of motion MM 



7](t) ± r?(r) + t] 3 {t) + 5 S(r) tj(t) = 



(3.10) 



— ±l + q 2 + V 2 (T)+gZ(T) 



(3.11) 



Two different cases correspond to the different signs in the evolution equations above. 

The negative sign is associated with tree level potentials that allow global 0(2N + 1) broken symmetric 
ground states, whereas the positive sign determines a potential with a symmetric minimum. As it will 
be discussed in detail below, the non-equilibrium dynamics in the broken symmetry case is described for 
early times by the process of spinodal decomposition and phase ordering and triggered by long-wavelength 
instabilities just as in a typical second order phase transition during a rapid quench through the critical 
temperature 

For positive sign, the physical situation that we want to describe is the case when the order parameter 
has an initial value corresponding to a large amplitude <p(t = 0) of order 0(m R /y/X R ) i.e. T]{0) — 0(1) [see 
eq. ( |3.8| )1. The subsequent non-equilibrium evolution of the order parameter is described in terms of large 
amplitude oscillations around the minimum of the potential. This situation would describe the dynamics 
after the phase transition when the order parameter has rolled down the potential hill and undergoes large 
amplitude oscillations near the minimum. In cosmology this situation also descri bes the period of reheating 
in chaotic scenarios (2jj2^|. As can be seen from the equation of motion ( 3.1l[) th e effective mass for the 
charged field modes oscillates in time leading to parametric amplification |24| , |25| , 136| . 

In this case the phenomenon is that of energy transfer from the 'zero mode' i.e. from the expectation value 
of the order parameter to the modes with non-zero wavevectors as a consequence of parametric amplification 
of quantum fluctuations. 
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Thus, the physics is very different between the two cases and the only feature in common is that either 
through the growth of long-wavelength fluctuations via spinodal instabilities or the growth of fluctuations via 
parametric amplification the ensuing non-equilibrium dynamics results in the production of a dense plasma 
of charged particles strongly out of equilibrium. 

The initial conditions on the order parameter (condensate) are chosen to be 

77(0)^770 , r)(0) = 0, (3.12) 
and the initial conditions on the mode functions are [see eqs.(^5|) and (|Q|)] 

<p q (Q) = -^= , <p q (0) = -iyftT q , (3.13) 
where the dimensionless frequencies Q q will be determined in each particular case below. 



A. Broken symmetry: spinodal instabilities 



Consider the case in which the system is undergoing a sudden phase transition out of equilibrium from 
an initial disordered state at large temperature very rapidly to almost zero temperature, i.e. a quenched 




For r > the equations of motion are those for 
with 77(0) <C 1 , 77(0) = 0. For simplicity, we shall 



phase transition with a vanishing order parameter 
a broken symmetry case with the (— ) sign in eq 
consider the case 77(0) = which entails t/(t) = 0. 

Furthermore, we s ee that for very weak coupling and early times i.e. when the back reaction from the 
term <?£(r) in ( 3.1 lj) can be neglected, there is a band of spinodally unstable wave- vectors < q < 1. The 



modes in this unstable band will grow exponentially initially. Because we are describing an initial condition 
corresponding to a sudden quench, we impose the initial condition that at the initial time the mode functions 
describe particles of the stable phase, i.e., we choose the initial frequencies for the modes in the unstable 
band to be given by ||3|$3|57j 



tt q = vV + 1 for q 2 < 1 
the short wavelength modes are not affected by the sudden quench and we choose 



1 for 



? 2 >1 



(3.14) 



(3.15) 



However, we emphasize that detailed numerical studies reveal that the dynamics is not very sensitive to the 
choice of the initial frequencies for weak coupling |p6[p7| . 

The important feature is that this initial state has non-perturbatively large energy density, of order 
\rn,Ft\ A I \r as compared to the broken symmetry vacuum state, for which |?7| = 1. 



As discussed in refs. |36 37 the ensuing dynamics is strongly out of equilibrium. The modes with wavevec- 
tors in the unstable band begin growing exponentially and their contribution to the self-consistent expectation 
value gS(r) causes it to grow exponentially. This instability is the hallmark of the process of phase separation 
and determines the emergence of correlated regions |34|,|3(||3lj : these are the familiar spinodal instabilities as- 
sociated with the process of phase separation and phase ordering. The contribution of these unstable modes 
to <?£(t) dominates the early time dynamics and when <?E(t) becomes of 0(1) and competes with the tree 
level term (—1) in the evolution equations for the mode functions (3.11) these instabilities shut-off through 
the backreaction. This defines a new dynamical time scale that determines the onset of full non-linear 
evolution and is estimated to be [t56|]37ri 



tnl = - In 



0(hx\lng\ 



(3.16) 



Thus, two different regimes emerge: 
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i) the early time regime for r < tnl hi which the back-reaction can be neglected and the evolution of the 
mode functions is essentially linear and dominated by the spinodally unstable wave-vectors for which the 
mode functions grow exponentially (linear instabilities). 

ii) the late time regime for r > t^l for which the effective mass squared A4 2 (t) = — 1 + <?£(t) + ?y 2 (r) 
tends to zero and the mode functions become effectively massless |57j . 

1. Early time regime 

For t < tjvl and weak coupling, the effects of the back- reaction can be neglected and the mode functions 
obey a linear equation of motion. Whereas the modes outside the spinodally unstable band oscillate and 
their amplitudes remain bound in time, those in the unstable band grow exponentially. For the case r] Q <C 1 



we can neglect at early times both the quantum fluctuations g£(r) and t] 2 (t) in the mode equations (3.11). 
The explicit solution is thus |37j 

ip q (r) = a q exp (ry/l - q 2 ^j + a* exp (-Ty/l - q 2 ^j (3-17) 



where the coefficient a q is determined by the initial conditions (3.13), i.e. 

^(0) = (l + g 2 )- 1/4 , ¥> 9 (G) = -<1 + <7 2 ) 1/4 , 

we find 



y/l-g 2 - yi + g 2 
2 V / 1^(1+9 2 ) 1/4 ' 



a i = „ r, z u ■ ( 3 - 18 ) 



A feature of the solution (3.17) with ( |3 .18 ) that will become important is that when the exponentially 



damped solution becomes negligible as compared to the exponentially growing one, the phase of the mode 
functions tp q (j) freezes, i.e. becomes constant in time and is a slowly varying function of q for long wave- 
lengths. 

2. Late time regime 

For times r > t nl the effective mass term M 2 {t) = — 1 + .gS(r) vanishes leading to the sum rule |3(|[37j 

gS(oo) = 1 (3.19) 

and the mode functions obey a massless wave equation. The asymptotic solutions are given by |3j]] 

tp q (r) = A q + Bq e-^ (3.20) 

where the coefficients A q , B q are both non- vanishing because the Wronskian is constant and determined by 
the initial conditions 

- 2* = W[p q , - <&,(r)^(r) - <p q (T)<p* q (r) = -2iq [\B q \ 2 - \A q \ 2 \ (3.21) 
leading to the important result 

\B q \ 2 -\A q \ 2 = I (3.22) 



Furthermore, the sum rule ( [3.19 ) is asymptotically dominated by the modes in the unstable band 
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fliS(r) T 3 / g 2 d<? [\A q \ 2 + \B q \ 2 ] + oscillating terms 
Jo 



(3.23) 



where the oscillating terms vanish as 1 jr. We conclude |37j that for the modes in the unstable band 

\A q \ 2 = 0(l/g) = \B q \ 2 , 0<g<l (3.24) 

determining that A q , B q are both of 0(1/-Jg) whereas A q , B q are of order one elsewhere. 

The following sum rules arise from a) the vanishing of the effective mass and b) conservation of energy 

13 

/ q 2 dq\A g \ 2 = 1 + 0(1) (3.25) 
Jo A 9 



f V dq \ Aq \ 2 = ^ + 0{l) 
Jo 8 ff 



(3.26) 



Furthermore, the small q behavior of A q and B q is given by |37| ], 

4r *=° -^[K + iqL + 0(q 2 )] , B q q =° ±[K- iqL + 0(q 2 )] (3.27) 

where ImLK = 1 according to the Wronskian condition |}7) . 

The non-zero coefficient K is determined by the linear growth in time of the mode <£a = o( r ) m this case 
with broken symmetry |j37f . For small coupling it is found numerically to be given by |37[ , 

K = K + — iK- where K± = — {l.l ... T 0.003 ... g + Dig 2 )} . (3.28) 

This asymptotic behavior for small momentum will prove to be important for a quantitative analysis of the 
magnetic mass. 

B. Unbroken symmetry: parametric amplification 

In th e unb roken symmetry case, corresponding to the choice of the plus sign in the equations of motion 
(|3.10D- (|3.1lD the frequencies fl q are chosen to be 



ft g = vV + l + J7 2 (0) (3.29) 

and the initial condition for the dimensionlcss order parameter is chosen to be 

77(0) = 770 = 0(1) ; r)(0) = (3.30) 

In this case the 'zero mode' (expectation value) 7?(r) oscillates around the minimum of the potential resulting 
in an oscillatory time dependent mass term for the modes ip q (r). 

1. Early time regime 

Neglecting the backreaction of the fluctuations, an oscillatory time dependent mass leads to parametric 
amplification of the mode functions which are Floquet solutions. These solutions are characterized by 
parametric instability bands. 

For weak coupling the early time behavior of 77 (r) an d the mode functions ^p q {r) can be found by neglecting 
the backreaction terms in the equations of motion ( [3.1C| )-( [3.11 ) for the unbroken symmetry case. The 
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equation for the zero mode with the initial conditions ( 3.30 ) has as solution a simple elliptic function p7[ . 
Inserting this elliptic function, the evolution equation for each mode fq{r) becomes a Lame equation that 
can be analytically solved in terms of Jacobi theta functions, the details are given in p7| . The important 
feature is that this Lame equation has only one band of parametric instability for real q. The unstable band 
corresponds to wavevectors |37| 



(3.31) 



The modes in the unstable band grow exponentially in time, whereas those in the stable region 770 /v2 < 
q < 00 oscillate in time with constant amplitude. 

The explicit solution with boundary conditions ( 3.13] ) for the mode functions in the unstable band is given 

by 



ip q (r) = a q U q (-r) + a* U q {r) 

with the Floquet solution U q (r) g iven in [ p7| . 

With the choice of frequencies ( |3.29 ) the coefficient a q is found to be given by 



(3.32) 



1 



2, AT 



2in q 



; W q = -2q 



l rt/2 + l + q* 
Til/2- q* 



(3.33) 



The Floquet solutions U q {r) are derived in detail in |31]] and depend on the initial condition rj through the 
nome (7(770). Since in this case (7(770) < e _7r = 0.0432139... for any initial condition 770 p7| , we can express 
g(?7o) by the excellent approximation 



1 (1 + 77 2 )V 4 - (1 + ^/2)^ 

™ 2(l + % 2 )V4 + (l + r? 2 /2)1 /4 



(3.34) 



with an error smaller than ~ 10 7 [j37| . In addition we can use the approximation 5(770) <C 1 and the Floquet 
solutions simplify in this limit to 



U q (-r)=e B ^ 



iin (nv q - V 1 + Va r ) 



sm ttv 



0(q) 



(3.35) 



with Floquet index 



B q = A^l + rf Q q(r] ) sin 2wv q + 0(f) , simrv q = , 1 - —q 2 + O(q) 

V % 



COS TTV n 



= ^q + 0(q) 



(3.36) 



Therefore the backreation gS(r) grows exponentially at early times because of the parametric instabilities. 
The exponential envelope of the backreation term is given by p7| 



9 e*' 



(3.37) 



where N and B can be found in |37j. When the backreaction competes with the tree level term, i.e. 
S'E(t) »1 + '7o/2 the full nonlinearities must be taken into account, this equality determines the non-linear 
time scale t^vl given by fl37fl 



TNL 



B 



N{l + r,l/2) 



gVB 



(3.38) 
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Detailed analytic and numerical studies in j37| reveal that most of the particle production occurs during the 
time interval r < tnl- 

For ttvl > t > 1 the modulus squared of the mode functions |<p g (r)| 2 is peaked at the value of q at which 
the Floquet index is maximum, this value is given by |}7j 

q* = \vo[l- q(m)} + O(f) (3.39) 



2. Late time regime 



The parametrically resonant band, < q < is shut-off by the non-linearities [the term gTi(r)] for 

times r > tnl. Two non-linear resonant bands appear in this regime. One near q — and the other just 
below q = r]o/V2- The width of these nonlinear resonances diminishes in time. We have for the non- linear 
resonant bands [BTl 



0<(? 2 <^I and ^ 
t 2 



K 2 . 2/^0 

- <q <y 



(3.40) 



(with Ki , K2 determined in ref . |37| ) and the phase space for these small resonant regions becomes increas- 
ingly smaller at late times. 

Asymptotically, the effective mass oscillates around the constant value |37| 



M 2 {oo) = 1 + gS(oo) = 1 + 



1o 



(3.41) 



and the mode functions can be written as 

V q {T)=A q {r) e^ + B q (r) e 



= vV+-M 2 (oo) 



(3.42) 



where the amplitudes A q (r) and B q (r) depend on the slow time scale t/t^vl for r > tnl and are defined 

by [p7l 



Mr) 



1 



Mr) = 2 e+M?T 



M T ) 



M T ) 



M T ) 



„• M T ) 



(3.43) 
(3.44) 



These amplitudes vary slowly in time and in particular ajg|Ag(r)| is identified with the number of asymptotic 
particles of mass A^ 2 (oo) |37| . 

For wavevectors inside the small bands of non-linear resonances ( 3.4C ), these amplitudes grow with a 
power law p?| , whereas the modes outside from these resonant regions oscillate with constant amplitude. 
The fact that the width of these resonances diminishes at longer times is a consequence of the non-linearities. 
A very important consequence is that asymptotically for all modes with q ^ 0, ri /y/2 



lim A q (r) = A q 



lim B q (r) = B q 



(3.45) 



where A q and B q are constants. 

Hence, the mode functions with q ^ , r/o/^/2 asymptotically behave as 

<p q (r) T = 1 A q e l ^ T + B q e"^ T 
Asymptotically the constancy of the Wronskian leads to 



(3.46) 
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\B q \ 2 -\A q \ 2 = - (3.47) 

Furthermore, for the modes with wavectors in the resonant band < q < rio/V^, 

\A g \ = o(^j=\B q \ (3.48) 

Just as in the broken symmetry case, there are two important sum rules as a result of the asymptotic value 
of the effective mass and of conservation of energy. In this case these read H^j 

rr)o/V2 -i /•»?o/\/2 i 

J dq q 2 \A q \ 2 = Tg Vl + (g°) , jf dq g 4 \A q \ 2 = — 4 + (g°) , (3.49) 



C. Formation of the plasma 



The main conceptual feature that emerges from the summary above is that in both situations, broken 
or unbroken symmetry, spinodal or parametric instabilities lead to profuse particle production. The parti- 
cles that are produced are charged scalars, these are produced in pairs of total zero momentum, and the 
distribution of produced particles is localized in the region of instabilities. In the case of broken symme- 
try the distribution is peaked in the region < q < 1 and in the case of unbroken symmetry in the region 
< q < r]o/y/2. In both cases the amplitude of the mode functions in these regions become IVgMI 2 = 0(\/g) 
i.e. non-perturbatively large. This amplitude is associated with the number of particles created |56]]37]] (see 
below) and therefore we conclude that during the period of spinodal or parametric instabilities < r < tjvl 
a dense plasma of charged particles is formed as a result of these instabilities. This plasma is neutral 
and is described by the distribution functions of the particles, which is proportional to \ip q (r)\ 2 fS6p^ ] (see 
below) and is clearly a non-equilibrium distribution in the sense that it cannot be described by a thermal 
distribution at some temperature. These distribution functions had been obtained numerically in p6| , p^7| . 
Figs. [I] and |^ display g\<p q (T~Ni J )\ 2 for the broken and unbroken symmetry cases respectively, it is clear that 
the square of the mode functions become of order 0(1/ g) at r « tml for wavevectors in the unstable bands. 

Furthermore the distribution of particles continues to evolve for r > tml and this evolution is more marked 
in the unbroken symmetry case. 

We now turn to the description of the electromagnetic properties of this non-equilibrium plasma of charged 
particles for which we need the non-equilibrium Green's functions of the charged scalar fields. 



D. Basic ingredients: real time non-equilibrium Green's functions 

The proper description of real time non-equilibrium evolution is in terms of the time evolution of an initial 
density matrix. A formulation in terms of a path integral along a complex contour in time allows to use the 
familiar tools of quantum field theory to study non-equilibrium phenomena. In this, the Schwinger-Keldysh 
or CTP (closed time path) formulation p3|]36| ], the essential ingredients are the non-equilibrium Green's 
functions. In particular there are four possibile two-point functions, denoted by indices (a, b) G {+, — } 
which correspond to the evolution along the forward and backward time branches. 

• Transverse photon propagators 

Since photons will be treated perturbatively, we need the bare photon Green's functions. Furthermore 
we will consider that the initial state is the photon vacuum. Therefore the relevant real time Green's 
functions for transverse photons are given by 

(4^,i)4 b j(zV)> = ~ i j-§^ s#(m,o e-*(*-*') , 
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where the explicit form of GfHk; t, t') is 



gt j + (k-,t,t')=v ij (k) [g>(t,t')Q(t-t')+g£(t,t')e(t'-t)] , 
^-(fc;i,t')=^(fc) [GZ(t, t')e(t' - t ) + g<(t,t')e(t-t')} , 
g±jk-,t,t') =v ij (k) g<{t,t>) ■ gm;t,t , ) = r ij (k) g>{t,t>) 



(3.50) 
(3.51) 
(3.52) 



and Vij(k) is the transverse projection operator, 

V ij {k) = 5 ij 



k^ kj 



(3.53) 



At tree level 



(3.54) 
(3.55) 



Scalar propagators 

The scalar propagators are truly non-equilibrium and can be written in the general form 

d 3 k 



(2tt) 



Gf(t,t') e -* k < s -^ 



where (a, b) S { + , — }, r, s = 1, . . . , N. With the field expansion give n by (3.2)-(3.3) in terms of the 
(dimensionful) mode functions f k (t) obeying the equations of motion (3^) we obtain in the N = oo 
limit, 



G++(t, = G>(t, t')0(t - t>) + G< (t, t')Q(t' t) , 
Gr{t, t') = G> (t, t')Q{t' -t) + G< (t, t')Q{t - if) , 
G+-(M') = G<(M') ; G^(t,t') = G>(t,t'), 

G > k (t,t')= l -f k {t)r k {t>), 



(3.56) 

(3.57) 
(3.58) 

(3.59) 



An important property of the mode functions in the asymptotic region will allow us to establish a cor- 
respondence between the non-equilibrium results to be obtained below and the more familiar equilibrium 
results. In both cases, broken or unbroken sym metry, after the non-linear time scale the mod e fun ction s 
become those of a free field theory [see eqs.( 3.2C ) and (|3.46|) 1. The Wronskian conditions ( 3.22 ) and ( 3.47 ) 



allows to write the modulus of the coefficients A q , B q in the form 



\B q \ 2 



[i+Ag 



\A q ? 



(3.60) 



with w q = q for broken symmetry or uj q = \/ q 2 + Al 2 (oo) for unbroken symmetry. 

Af q describes the distribution of asymptotic charged scalar particles created during the rapid non- 
equilibrium stages of spinodal decomposition or parametric instabilities |^6| , ^7| . These are non-equilibrium 
distribution functions, a result of profuse particle production during the stage of spinodal instabilities or 
parametric amplification. The number of charged scalars produced during these stages is non-perturbatively 
large, since for the wave vectors in the unstable bands J\f q is of order 1/g. 
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The asymptotic behavior of the functions G k (t, t') when both time arguments are in the asymptotic region, 
much larger than the non-linear time scale can be written in the illuminating form 



G>{t,t') T = — (1+JVfc) e-*"^*-*) +JV k e^C"*') +iR e A k B* k e lUk{ - t+t,) 
2uj k L J L 

G<{t,t') = G>{t',t) (3.61) 

where the mixing terms proportional to e ±luJkf - t+t ) are a signal that the non-equilibrium behavior remains in 
the asymptotic region. In most circumstances these rapidly varying oscillatory terms lead to contributions 
that vanish very rapidly by dephasing. 

The first two terms of the function ( |3.61| ), that depend on the difference of the time arguments can be 
compared to that of a free field theory in equilibrium 



k.equil 



M = ^{ [1 + Hk] e-^-^+rike^-t')} ; G< equil {t,H) = G>(f,t) (3.62) 



where n k is the thermal distribution function. Thus we see that the part of the asymptotic non-equilibrium 
Green's functions that depend on the difference of the time arguments has the form of the free field equilibrium 
Green's functions but in terms of the non- equilibrium distribution functions J\f k . This formal similarity will 
allow us to compare the non-equilibrium results in the asymptotic regime to those more familiar in equilibrium 
field theory and to interpret the different processes in the medium. 

It is convenient to summarize the main features that will be responsible for the phenomena studied below. 

• In either case, spinodal instabilities in the case of the potential allowing broken symmetry states or 
parametric instabilities associated with the oscillatory evolution of the order parameter in the case of 
unbroken symmetry, there are strong fluctuations that lead to non-perturbative particle production 
of the order (D(l/g). At times larger than the non-linear time, when the effects of the backreaction 
become of the same order as the tree level terms, the occupation number of modes in the unstable 
band is non-perturbatively large [0(1/ g)]. The state of the system can be best characterized as a 
non- equilibrium dense plasma. The distribution function of the created particles is not an equilibrium 
one and has a finite limit for infinite time. 

• As a result of non-perturbative particle production, the Green's functions of the scalar fields, deter- 
mined by the mode functions f k (t), are those of a plasma strongly out of equilibrium and will provide 
a non-perturbative contribution to the photon polarization. 

• To leading order in the large N limit and to lowest order in a = e 2 /47r, the photon polarization is 
given by the diagrams shown in fig. ^a,b. 

The loop is in terms of the full scalar propagator in the leading order in the large N limit, which 
receives contributions from the mean-field and background as depicted in fig. ^p. 

We now have all of the ingredients to study the electromagnetic signatures of these non-perturbative 
phenomena to leading order in the large N limit and to lowest order in a — e 2 /An. 

We emphasize again that these phenomena have nothing to do with the ordinary Higgs mechanism. In this 
model the global gauge symmetry is not spontaneously broken by the initial state even when the potential 
for the scalar fields allows for broken symmetry and remains unbroken throughout the dynamical evolution. 



IV. PHOTON PRODUCTION VIA SPINODAL AND PARAMETRIC INSTABILITIES 

We study the production of photons both via the spinodal instabilities associated with the process of phase 
ordering (in the case of broken symmetry potentials) and via the parametric instabilities associated with the 
non-equilibrium evolution of the order parameter around the symmetric minimum (in the case of unbroken 
symmetry potentials). As described in the previous section, we carry out this study to leading order in the 
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large N limit and to lowest order in the electromagnetic coupling. This is similar to the formulation in 
]4l| , |f2|| for the rate of photon production in the QGP to all orders in a s and to lowest order in a em . 

However, our approach differs fundamentally from the usual approach in the literature p|-pl5||, which 
relies on the computation of the photoproduction rate from processes that satisfy energy conservation, i.e. 
on shell. This results from the use of Fermi's Golden Rule in the computation of a transition probability 
from a state in the far past to another state in the far future. 

Instead our computation relics on obtaining the integrated photon number at finite time t from the time 
evolution of an initial state at tg. Clearly this approach is more appropriate in out-of-equilibrium situations 
where transient, time-dependent phenomena are relevant. 

Non-equilibrium time dependent transient phenomena cannot be captured by the usual rate calculation 
based on Fermi's Golden Rule, since such calculation will obtain the number of produced photons divided 
by the total time t in the limit when t — > oo. This definition is insensitive to the non-energy conserving 
processes, which are subleading in the long time limit but could dominate at finite time, and could potentially 
lead to grossly disparate estimates of the total number of photons produced in a situation in which a plasma 
has a finite lifetime as is the case in heavy ion collisions. 

To lowest order in a and leading order in the large N the leading process giving rise to photoproduction 
is the off-shell production of a pair of charged pions and one-photon from the initial vacuum strongly out 
equilibrium. Thus, we consider the transition amplitude for the process |0 >— > |7r + 7r~7 > to order e, more 
precisely the amplitude for the interaction to create a pair of scalars with momentum q + k and q respectively 
and a photon of momentum k and polarization A. The initial state |0 > at time to is the Fock vacuum for 
pions and photons but its evolution is non-trivial because it is not an eigenstate of the Hamiltonian, nor 
is it perturbatively close to an eigenstate. In the case of spinodal instabilities this state is unstable, and it 
decays via the production of pions and photons. In the case of parametric instabilities, this state involves a 
dynamical expectation value for the a field with non-perturbatively large amplitude (i.e, r)o ~ 1). 

The lowest order contribution to this amplitude in the electromagnetic coupling is given by 



t 

3„ 



to 



A,fe,A(*) =< 7T 7T l\ i I dtidrxJ(ti,x)-A T (ti,x)\0> (4.1) 
where J(ti,x) is the electromagnetic current 

JV 

J(tx, x) = 4= E (^ V ^ - V ^ $ ( 42 ) 

V r— 1 



If the (transverse) photon field is expanded in terms of creation and annihilation operators of Fock quanta 
associated with the vacuum at the initial time 

l T (t,x) = E / e~ ikt+i ^ a x (k) + (e x (k) e -«+ ifc ^)* a{(k)] (4.3) 



X=l,2 



and the scalar fields are expanded as in eqs.(p.2[)-(p.3|), we find the amplitude to be given by 



A q , k , X (t) = ^== J t dh e x (k) ■ q T e^ r q {h) r mi] (h). (4.4) 



ikt\ £ * 

uli CAW ' e 

to 

Squaring the amplitude, summing over q and r and A and using 

2 



E4(£) e{(k)=r^(k) 



x=i 



we finally obtain that the total number of photons of momentum k produced at time t per unit volume from 
the initial vacuum state at time to is given by 
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N ph (k,t) = (2tt)' 



<i 6 iV(£) 
d 3 ir <i 3 fc 



e 

2A 7 



rf 3 g 



9 2 (1 



cos 



/ 9 (^i)/| g - + fc|(ii)e- Jfctl dt x 



(4.5) 



where 6 is the angle between q and k, q- k = q k cos 0. The same formula can be obtained as a particular case 
of the generalized kinetic equation for the photon distribution function obtained in Appendix A. We refer 
the reader to this Appendix for a more detailed discussion of the kinetic equation and its regime of validity. 



We point out that if the mode functions f q {t) are replaced with the usual exponentials exp(- 



OA 



and the limits to — > — oo and t — > oo are taken, the familiar energy-conservation Dirac delta function is 
recovered and therefore the process is kinematically forbidden in the vacuum. Furthermore, the discussion 
in the previous section highlighted that during the stage of spinodal instabilities or parametric amplification, 
the mode functions in the unstable bands grow exponentially. Hence the modes in the unstable band will lead 
to an explosive production of photons during these early stages. Clearly the maximum production of photons 
will occur in the region of soft momenta, with the wavevector k in the unstable bands. In this manner the 
scalar mode functions with wavevectors q and q + k will be in the unstable bands leading to four powers 
of the exponential growth factor. Thus we will focus on the production of soft photons studying the case 
of broken symmetry (spinodal instabilities) and unbroken symmetry (parametric instabilities) separately. 
Having recognized the emergence of a dynamical time scale 



TNL 



mil 

,9 



[see eqs.(3.16) and ( 3.38| ) for more detailed expressions] that separates the linear from the non- linear behavior, 
we analyze both regions r < tnl and r > tnl separately. 



A. Photon production via spinodal instabilities: 



1. r < t nl 

The number of produced photons of wavelength k per unit volume at time t is given by eq.([ll]). Obvi- 
ously the integrals in this expression can be computed numerically Q since the mode functions are known 
numerically with high precision ^T\ . However, the summary of properties of mode functions for r < tnl and 
t > tjvl provided in the previous section allows us to furnish an analytical reliable estimate for the photon 
production. During the early, linear stages, we can insert the expression for the mode functions given by 
eqs. (|3~T7|)-(|3~18|). Furthermore, we focus on small k so that q and q + k are in the spinodally unstable bands 
and keep only the exponentially increasing terms which dominate the integral at intermediate times. The 
time integral can now be performed and we find (using dimcnsionlcss units) 



N ph (k,T) = — 



2k J (2tt) 3 



exp 



q 2 (1 - cos 2 6>) \a q aq+kl 2 



2 T 



V^~7+^Ji-(q-+kr 



V^~7+^/i-(q + ky 



(4.6) 



+ k 2 



Furthermore, the dominant contribution to the q— integral arises from the small q region justifying the 
non-relativistic approximation q 1, q + k 1. Hence N p h(k, r) becomes 



i_£Lr(4-fc 2 ) f°° dq 
16 2k J (2tt)2 

For k 2 r 3> 1 we can use the approximation 



dx q 4 (1 - x 2 ) e - 2 ^ 2 +^) 



(4.7) 



dx (1 - a; 2 ) e- Ax = 2e 



T2 +0(l/A*) 
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to perform the angular integration. Notice that the dominant region corresponds to x = — 1. That is q and 
k in opposite directions. Physically, this corresponds to two charged scalars with parallel momenta q and 
q — k emitting a collinear photon with momentum — k (see fig. 
In this regime the photon spectrum becomes 

1 ,= 2 f°° i 

" = ss*s i * d^r- v (4 - 8) 

where the factor l/(/cr) 2 arises from the angular integration. 

Now it is possible to compute the momentum integral via a saddle point approximation. Using the saddle 
point q = q = k/2 we obtain 

N v^r) = ^ e ^/2 + o (fc *)] + c(1/r) (49) 
where the proportionality factor iVi is given by 

V2 



m = 



2048 7T 3 / 2 ' 

We see that for r < tnl the number of produced photons grows exponentially with time. The production is 
mostly abundant for soft photons k <C 1. However, the derivation of (^^) only holds in the region in which 
the saddle point expansion is reliable, i.e. for k 2 r ^> 1. 

The k — > limit can be studied directly. In such case the angular integration in eq.(4.7) is straightforward 
and the momentum integration can be done using the result 



/>oo 

/ dq <7 4 exp(-2T<? 2 ) 
Jo 



'2tt 



64 r 5 / 2 



leading to the same result as eq.(4.9), 

We thus find an exponentially growing number of emitted photons (as ~ e 4r ) for r < tnl- Since 
e ir NL ^ we gee the total number of emitted photons at r k tml is of the order 

e 2 

N ph (k,t NL ) ~ — ^ (4.10) 
K g 

and is predominantly peaked at very low momentum as a consequence of the fact that the long-wavelength 
fluctuations are growing exponentially as a consequence of the spinodal instability. The power spectrum for 
the electric and magnetic fields produced during the stage of spinodal growth of fluctuations is 

(\E(k,r)\ 2 ) « (\B(k,r)\ 2 )^kN ph (k,r) . (4.11) 

Two important results can be inferred for the generation of electric and magnetic fields 

• At the spinodal time scale r « t^l the power spectrum is localized at small momenta and with 
amplitude ~ a/g 2 . 

• Taking the spatial Fourier transform at a fixed given time we can obtain the correlation lengt h o f 
the generated electric and magnetic fields. A straightforward calculation for r < tml using eq. (4.9) 
reveals that 

(E(r,r) ■ E(6,t)) ~ (B(r,r) ■ B(6,t)) ~ e"i^ ; £(r) ~ (4.12) 

The dynamical (dimensionful) correlation length £Jt) ~ y/r is the same as that for the scalar fields 
before the onset of the full non-linear regime |3^,^]. Therefore, at early and intermediate times the 
generated electric and magnetic fields track the domain formation process of the scalar fields and reach 
an amplitude ~ a/g 2 at time scales r tjvl over length scales ~ |to.r| -1 [k^l/g)] 1 ^ 2 . 
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2. T > TNL 



We now split the time integral in eq. Q4.5| ) into two pieces, one from up to tjvl and a second one from 
t~nl up to r. In the first region we use the exponentially growing modes as in the evaluation above, and in 
the second region we use the asymptotic form of the mode functions given by eq.( |3.46 ). The time integral 
in this second region can now be performed explicitly and we find 



N ph (k,T) 



8tt 2 k 



q 4 dq 



dx (1 - x 2 ) 



A 1 B \q+k\ 



B 1 B \q + k\ 



dri v q (n) <p^(n) e ^ + A q A w+K{ 

e i(q-\q+k\-k)r _ J.{q-\q+k\-k)T N t 

q — \q + k\ — k 

p -i(q+W+k\+k)r _ -i(q+\q+k\+k)r NL 



j(q+\q+k\-k)r _ i(q+\q+k\-k)r NL 



q+\q + k\+k 



B 1 A \q+k\ 



q+ \q + k\ - k 

e -i{q-\q+k\+k)T _ e ~i(q~ W+R\+k)T N1 



q-\q + k\+k 



(4.13) 



where \q + k\ = \J q 2 + k 2 + 2 k q x. The momentum integration is restricted to the region of the spin- 
odally unstable band since only in this region the modes acquire non-perturbatively large amplitudes. The 
integration over q > 1 only provides perturbative O(g) corrections. 

The contribution of the asymptotic region r 3> tnl m eq.(4.13) displays potentially resonant denominators. 
As long as the time argument r remains finite the integral is finite, but in the limit of large t 3> tnl the 
resonant denominators can lead to secular divergences. In the long time limit we can separate the terms 
that lead to potential secular divergences from those that remain finite at all times. Close inspection of 
eq.(4.13) shows that asymptotically for large time the square modulus of the second, third and fourth terms 
yield potential secular divergences. The square modulus of the last term is always bound in time and 
oscillates since the denominator never vanishes. In addition, the cross terms either have finite limits or are 
subdominant for t — * oo. The square modulus of the first term is N p h(k,tffL) given by eq. (4.10). In order 
to recognize the different contributions a nd to establish a relationship with the equilibrium case it proves 
useful to use the definitions given in eq. (3.60). We find the following explicit expression for the dominant 
contributions asymptotically at late times, 



N p h(k, 



t) T = l N ph (k,T NL ) 



At: 2 k 



q 4 dq 



(q+ W+k\ 



q \q 

(r- 



+ i 



dx (1 - x 2 ) x 



TNLi 



q+\q + k\-k 



1 — cos 



q- \q + k\ - k) (r- t nl ) 



q — \q + k\ — k 



+ [l+M q ] M w+%] 



1 — cos 



q - \q + k\ + k) (r - t nl ) 



q - \q + k\ + k 



>[l + O(g)]+O{T ) 



(4.14) 



The first term, containing the factor jV g A/j- + ^|, corresponds to the process tt + tt~ — * 7, i.e. massless 
charged scalar annihilation into a photon, the second and third terms (which are equivalent upon re-labelling 
q — * — q— k) correspond to bremsstrahlung contributions in the medium, tx^ — > + 7. 



The following relation 



1-x' 



1 



Aq 2 k 2 



q + \q + k\ 



q-\q + k\ + k) (q-\q + k\-k) [q+\q + k\+k 
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ensures that there are only simple poles in the inte grand of cq .(4.14). 

Asymptotically for long time the integrals in eq.( 4.13 )-( 4.14 ) have the typical structure 



^(1 

o y 



cosyt) p(y) p(0) log [fie' t]+ [ — [p(y) 

Jo V 



.p(O)0Ou-y)] + 0(~ 



(4.15) 



where p(y) is a continuous function, uan arbitrary scale and and 7 = 0.5772157 ... is the Euler-Mascheroni 
constant. Notice that the expression ( 4.15 ) does not depend on the scale /i, as can be easily seen by computing 
its derivative with respect to /1. 

Therefore the simple poles arising from the collinear singularities translate in logarithmic secular terms 
appearing for late times acc ordin g to eq .(4.15). 

The denominators in eqs.( 4.13 )-( 4.14 ) vanish leading to collinear singularities, i.e. kinematical configura- 
tions where the pho ton a nd a charged particle have parallel or antiparallcl momentum. More precisely, the 
denominators in eq.(4.14) vanish at the following points: 



\q + k\ = k — q, q — k and q + k , 

corresponding to cos 9 = x = — 1, —1 and +1, respectively. 

It is convenient to perform the angular integration using the variable £ = \q + k\ with dx = Since 
the most relevant contribution arises from the region of momenta inside the spinodally unstable band with 
J\f q = 0(1/ g) ^S> 1 the angular integration simplifies and we find 



N ph {k,T) 



t>1 



2tT fc 3 



log^i' 



qdq Af q [\q - k\ A^_ fe , + (q + k) N q+h ] [1 + (<?)] + O (r°) (4.16) 



Here \q — k\ stands for the absolute value of the difference between the numbers q and k. If we restore 
dimensions and we recall that M q is of order 0(1/ ' g) for < q < 1, we find that the logarithmic term has a 
coefficient ~ e 2 |m^| 3 /[g 2 fc 3 ]. This remark will become useful when we compare later to a similar logarithmic 
behavior in the case where the scalars are in thermal equilibrium (sec. VIII). 

These logarithmic infrared divergences lead to logarithmic secular terms much in the same manner as in ref. 
[ prOf and indicate an obvious breakdown of the perturbative expansion. They must be resummed to obtain 
consistently the real time evolution of the photon distribution function. The dynamical renormalization 
group program introduced in ref. p6[ provides a consistent framework to study this resummation. 

A similar logarithmic behavior of the occupation number has been found in a kinetic description near 
equilibrium in the hard thermal loop approximation pgj ] . 

Furthermore, we note that the evolution equation for the photon distribution function under consideration 
has neglected the build-up of population of photons, and therefore has neglected the inverse processes, such 
as charged-scalar production from photons and inverse bremsstrahlung. These processes can be incorporated 
by considering the full kinetic equation described in Appendix A. Hence a consistent program to establish 
the production of photons beyond the linear regime must i) include the inverse processes in the kinetic 
description and ii) provide a consistent resummation of the secular terms. We postpone the study of photon 
production in the asymptotic regime including these non-linear effects to a forthcoming article. 



B. Photon production via parametric amplification 



We now study the process of photon production during the stage of oscillation of the order parameter 
around the minimum of th e tree level potential in the unbroken symmetry case. This case co rresp onds to the 
evolution equations (3. 10)- (3. 11) with the plus sign and with the initial conditions fl3.29 )-( 3.30| ). We begin 
by studying the early time regime. 
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1. T < T N L 



The dominant contribution to the production of photons again arises from the exponentially growing terms 



in the parametrically unstabl e ban d. Hence we keep only the exponentially growing Floquet solution (3.35) 



with Floquet index given by (3.36) 



In order to perform the time integration we focus on the exponentially increasing terms and neglect the 
oscillatory contributions in the product 



Sin TTV. 



(^TTV q - yjl + rjl t \ sin (nv q+k - <J 1 + rj 2 r^j = sin 2 (irv q - <J 1 + rj 2 t \ + 0{k) 



oscillatory terms. 



In keeping only the exponentially growing contribution and neglecting the oscillatory parts we evaluate the 
envelope of the number of photons averaging over the fast oscillations. 



With these considerations we now have to evaluate the following integral according to eq.(4.5), for t»1, 
but t < tnl 



N ph (k,t) 



8k (2tt) 2 J 



dq / dx (1 — x 2 ) q 4 



Sill TTV q Sill TTVq+k 



exp[(B q + B q+k )T] 



(Bq + B, 



q+k) 



k 2 ' 



(4.17) 



As noted in |57J the Floquet index is maximum at q — r/o/2 and this is the dominant region in the q— integral. 
The fact that during the stage of parametric resonance the integral is dominated by a region of non- vanishing 
g is a striking contrast with the broken symmetry case and a consequence of the structure of the parametric 
resonance. As before, the strategy is to evaluate the integral for large times by the saddle point method. 
For q near ?yo/2 and for small k the saddle point is given by 

q = m /2 + x 2 k 2 +0(k 3 ), 

therefore the q— integral in the saddle-point approximation yields the result 



V 4 o 



256tt 2 k t 1 / 2 



exp 



r 16J1 



Vo 9 



o 



(4.18) 



For » 1 the angular integral (over x) is dominated by the region near x = and can be evaluated 
by using another saddle point expansion. In this limit the photon production process is dominated by the 
emission of photons at right angles with the direction of the scalar with momentum q. Physically this 
corresponds to two charged scalars with momenta q and q + k emitting a photon with momentum k with 
q-k = (sec fig. ||). This is another difference with the broken symmetry case wherein the production of 
low momentum photons was dominated by collinear emission. 

In this limit k 2 r 3> 1 the saddle point approximation to the angular integral yields the final result for the 
photon distribution function 



N ph {k,t) 



N 2 (m) 

k 2 T 



e 4T[B( V0 )+O(k 3 )] + [ - 



k 2 T > 1 



(4.19) 



where the coefficient B(j)q) in the exponential is given by 



B(vo) 



i + Vo q(vo) 



with q(rio) the nome given by eq. ( 3.34 ), and the factor N2(rjo) is given by 



32768^ (5 77 2 + 4)(3 7? 2 +4) 2 §3( % ) 



23 



and we note that an additional power r 1 I 2 in cq.(4.18|) arose from the angular saddle point integration. 



Wc find that there is a strong dependence on the initial condition of the order parameter tj(t), i.e. on 
rjo, which determines completely the energy density in the initial state. This is consistent with the strong 
dependence on the initial conditions of the mode functions that determine the evolution of the scalar fields 

In particular, we obtain for large 7]q 

2 

N p h (t) m = 1 * J* e° 2 "° T ; d = 373.83 ... , C 2 = 0.69142 . . . , k 2 r > 1 . (4.20) 
Ci k t 

Furthermore we also point out that in the region k 2 r 3> 1 there is an enhancement in the photon spectra at 
small momenta as compared to the broken symmetry case. This is a consequence of the photon emission at 
right angles (x = 0) in contrast with the collinear emission (x = ±1) for the broken symmetry case. 

In the range k 2 r < 1 the saddle point evaluation of the angular integral is not reliable, however in the 
very small k limit the angular integration can be done directly. We find 



N ph {k,t) k ^^01e^^^ + O(^ , 
where the proportionality factor N 3 (r] ) takes the value 

1 V2% 5 (l + % 2 ) 3/4 



k 2 r < 1 , (4.21) 



N 3 M 



~ 6144tt 3 /2 ( 5 v l + 4)(3 7? 2 + 4)2 q 5/2^ ■ 
In particular, for large 770 we obtain 

2 i/ 2 

N ph (r) "°= 1 * J*" e c ^" T ; C[ =422.60... , C 2 = 0.69142 ... , k 2 r « l . (4.22) 
This analysis reveals that the soft photon spectrum diverges as 1/fc and not as 1/k 2 at k — > . That 



guarantees the electromagnetic energy density (4.11) is infrared finite 



2. T > T NL 



For times r > tnl we use the asymptotic form of the mode functions given by eq. ( [3.42 ), we insert eq.( 3.42| ) 



in the expression (4J5) and we split the time integral into two domains < t\ < t^l and tjvl < n < r. The 



integral from tjvl < Ti < r is perf ormed explicitly with these asymptotic mode functions thus obtaining 



an expression analogous to eq.( 4.13 ). In this case, however, the upper limit of the momentum integration 
is Imax — ^0/^/2 i.e. the upper limit of the resonant band which gives the dominant contribution 0(1/ g 2 ) 
The integration over momenta q > q m ax gives a correction perturbative in g. We obtain , 



WV2 r+l 



I) 



TNL 







d n ^(n) y k - + S|(Ti) e lfcT1 + 

g i(w,-w |(r+£| -*:)T _ e i(u t -u w+nr k)T NL 



1 |g+fc| 

W 9 ^If+fcl + K 



24 



B 1 B \q+k\ 



^1 + + k 



[1+0(3)] 



(4.23) 



We focus on studying the small k behaviour < k -C 1 which can be obtained with the approximation 

k q x 



O (k 2 ) 



With this approximation the denominators in eq.(4.23) become, 



UJ q- UJ \<f+k\- k -- k l 1 



q x 



u W+k\ + k 



We rem ark that since .M 2 (oo) is non-zero, these denominators never vanish. Therefore the integrals in 
eq.( 4.23| ) do not generate secular terms and they have a finite limit for t — > 00. For as ympt otically long time 
and small k, the two denominators linear in k and their cross-product dominate eq.( 4.23| ). Isolating these 
dominant contributions we find 



N P h(k, 00) 



fe^O 



47T 2 k 3 



q 4 dq 



M q {l+M q ) F(q 2 )[l + 0(g)} 



(4.24) 



where F(q 2 ) is the regular function 



F(q 2 ) = J + dx (1 



2 3 + (g X/^) 2 = 2_Wg 

X ) [1 - (g i/w,)2]a q 



+ 3A4 2 (oo) 



ArgTh— - 



3w„ 



(4.25) 



With the identifications given by eq. (3.60) we recognize that the dominant contribution in the asymptotic 
regime to soft photon production arises from bremsstrahlung of massive charged scalars in the medium. 

A noteworthy feature is that the soft photon spectrum is strongly enhanced for small k since N p h(k, 00) 
grows as k~ 3 for k — * 0, this behavior must be compared to the d istribution at early time t < tnl where 
we had previously found that N p h(k,r) oc fc _1 for k — > [eq.( 4.21 )]. 

Thus in both cases, broken and unbroken symmetry, we find that the asymptotic non-equilibrium photon 
spectrum behaves for long wavelengths as f/fc 3 for k — > 0. This behaviour signals an IR divergence which 
may require a resummation of higher order terms in a. This is beyond the scope of this study. 

It will be found in section VIII that for charged particles in equilibrium, the photon spectrum has very 
similar features. Therefore, the total photon number N p h : TOT{T~) = J d 3 k N p h(k,T) is logarithmically 
divergent at small k. Nevertheless the total energy dissipated in photons, 



E p h,TOT(r) = J d 3 k k N ph (k,- 



is finite at finite times. As mentioned above, for late time in the broken phase, a resummation in a is needed 
to assess more reliably the photon distribution. 



V. PHOTOPRODUCTION FROM CHARGED SCALARS IN THERMAL EQUILIBRIUM 

We compute here the photoproduction process to leading order in e 2 from charged scalars in thermal 
equilibrium to compare it with the non-equilibrium case studied in sec. IV. However just as in the non- 
equilibrium case, we study the production of photons as an initial value problem, i.e, an initial state is 
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evolved in time and the number of photons produced during a finite time scale is computed. We emphasize 
again that this calculation is fundamentally different from the usual formulation of the rate obtained by 
assuming the validity of Fermi's Golden Rule and energy conservation. 

We shall find that there are some striking similarities between the two cases by identifying the high 
temperature limit T/\mn\ ^> 1 of the equilibrium case with the small coupling limit g -C 1 of the non- 
equilibrium situation. In both cases the plasma has a very large particle density. 

We highlight the most relevant aspects of the result before we engage in the technical details so that the 
reader will recognize the relevant points of the calculation. 

• Both from charged scalars in and out of equilibrium the photon production is strongly enhanced in 
the infrared since N p h(k,t) increases as 1/k 3 whereas for early times N p h(k,t) grows as 1/k. 

• In the broken symmetry case both in and out of equilibrium the number of produced photons increases 
at late times logarithmically in time due to collinear divergences. The physical processes that lead to 
photon production can be identified with collinear pair-annihilation and bremsstrahlung of pions in 
the medium. 

• The distribution of produced photons approaches a stationary value as t — > oo in the unbroken case 
both in and out of equilibrium with a distribution N p h(k, oo) ~ 1/k 3 . The relevant physical process is 
off-shell bremsstrahlung n — > tt + 7. 

Consider that at the initial time to there is some given distribution of photons N k (to) and charged scalars 
n p . The kinetic description provided in Appendix A leads to the following expression for the change in the 
photon distribution when the Green's functions of all fields are the form of the equilibrium ones given by eq. 
(3.62) but in terms of n p and Nk(to) [ p| 



N k (t)-- 



16T7 3 k 



d 3 C 



q 2 {1~X 2 ) / dt'{ 



cos[( Wg + u w+ ^ +k)(t- t')} [[1 + JV fc (t )](l + n,)(l + n w+ ^) - N k (t ) n q 
cos[(o; g + -k){t- t')] [[1 + N k (t )] n q n |ff£| - N k (t ) (1 + n q )(l + n |?+jf| ) 

cos[(w, - w |?f£| +k)(t~ t')} [[1 + JV fc (t )](l + n q ) n m%l - N k {t ) n q (1 + 
cos[(w g - -k)(t- t')] [1 + N k (t )]{l + n [f+ g|) n q - N k (t ) n |ffjE| (l + n q ) } 



(5.1) 



The different contributions in the above expression have a simple and obvious interpretation in terms of 
gain- loss processes |2q|. 



A. Photoproduction at first order in a 



In order to compare to the non-equilibrium situation described above, we will set the initial photon 
distribution to zero, i.e. N k (to) = 0, and we will also neglect the change in the photon population (this is 
also the case for the rate equation obtained by (i^^]). Integrating in time we obtain the expression 

N ph (k, t) = ^ J q 4 dq f dx (1 ~ [Ax +A 2 + A 3 + A 4 ] (q, k, x) (5.2) 

where 



J \q+k\ 



l-cos[ai(t-t )] 
Ai{q, k, x) = n q n w+%] — 2 , a x = uj q + 
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. , , s ri i 1 — cos[a 2 (t — to)] , 
A 2 (g,fc,x)=n, [l + n |f+jg| ] -2 , a 2 = u q - u,^, - ft . 

< JL 2 

a , i \ r-i i I - cos[a 3 (t — t )l , 
i43(g,*,aj) = [l + n,]n lstfjE | -g , a 3 = w, - + ft . 

< /■ i \ ri ,r- i 1 — COs[a 4 (t — t )] , 

M{q, k,x) = [1 + n q \ [1 +n^ + ^\ g ' "4 = w 9 + + ft . 

From this explicit expression one can easily see that in the zero temperature limit there is no photoproduction 
up to order e 2 . In fact, in the vacuum, only the term proportional to A 4 (q, k,x) and corresponding to the 
virtual process |0 >^ \tt + tt~j > remains but its contribution vanishes as 1/t in the long time limit, since 
the energy conservation condition 

a 4 (q, k,x) = uj q + + k = 

cannot be satisfied for positive non-zero u> q) u>^ + ^,k. This observation highlights that photon production 
will be completely determined by the plasma of charged scalars both in and out of equilibrium. We study 
in detail both cases separately. 



1. Broken symmetry phase 

In this case we study the spectrum of photons escaping from a thermal bath of massless scalars (Goldstone's 
bosons) with energy u q = q. The analysis is very similar to that performed in the non-equilibrium case and 
hinges upon extracting the secular terms in the asymptotic limit \i t 3> 1 . These arise from different kind of 
on-shell processes: 

1. the term Ai(q,k,x) corresponds to the annihilation tt + tt~ — > 7 in which a hard photon (ft > q) is 
emitted in the opposite direction of the initial pion {x = —1); 

2. the term ^2(9, ft, x) corresponds to the bremsstrahlung 7r — > 7r + 7 in which a soft photon (ft < q) is 
emitted in the opposite direction (x = —1); 

3. the term A%(q, k,x) corresponds to the bremsstrahlung tt — > ir + 7 in which the photon in emitted in 
the same direction of the pion (x = 1). 



Using eq.( |5.2D with uj q = q we recognize that the secular terms are of the same type as those of eq.( 4.15| ) 
and lead to a logarithmic divergence log /1 1 with /1 an infrared cutoff. After a detailed analysis similar to 
that carried out in the non-equilibrium case, we obtain 

t 1 e 2 f°° 
N p h(k 1 t) tlt = log fit J dqqn q [n lq _ kl \q-k\ + n q+k {q + k)] (5.3) 



which is remarkably similar to eq.(4.16) upon the replacement for the occupation numbers. For a thermal 



distribution of charged scalars the momentum integral is finite and for T 3> ft we find 



,2 rp3 



Nph (k,t) T ^ k j ^ log fit . (5.4) 

The high temperature limit of eq.(5.3) can be compared to the result out of equilibrium [cq.( 4.16] )] by 
identifying (T/m) 3 in the thermal case with 1/g 2 in the non-equilibrium case. In other words, m g~ z/ sets 
the scale of an 'effective temperature' to allow a qualitative comparison between the asymptotic description 
of photon production from charged particles with a thermal distribution and from a non-equilibrium plasma. 
However we emphasize that the non-equilibrium distribution is far from thermal and such a comparison 
only reflects a qualitative description. Furthermore, it becomes clear that the logarithmic secular term signals 
a breakdown of the perturbative kinetic equation and a resummation and inclusion of inverse processes will 
be required to study the long time limit. 
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2. Unbroken phase 



Also in this case the analysis is similar to the out of equilibrium computation: the final result is finite as 
t — > oo since there are no secular terms and we can simply neglect the oscillatory pieces. This is due to the 
presence of a non-zero mass for the scalars: as a consequence u> q = \/ q 2 + m 2 > q and the denominators in 
cq.( |5.2| ) never vanish (there are no collinear divergences). However for small k the two denominators linear 
in k 



a 2 {q, k, x) ~ -k(l + q x/ojg), a 3 (q, k, x) ~ k(l - q x/u q ), 



dominate and the formula simplifies as follows as t — > oo : 



N ph (k,oo) 



fe^O 



4-7T 2 k 3 



q dq 



+ 1 



f 



dx (1 — X 2 i \ — ■ ■ — - — ■ ; 

1 ^l^q+k [ V 1 + 1 X / U} q) (1-3 X /Uq) 



(5.5) 



From this expression one extracts a clear physical interpretation of the photoproduction process as generated 
by the off-shell bremsstrahlung of charged scalars in the medium. To give an estimation of N p h(k, oo) in the 
small k and high density limits we rewrite the previous formula as 



N ph (k, oo) 



An 2 k 3 



qUq F ^ 



n g (l+n q )[l + 0{k)] 



(5.6) 



where F eq (q 2 ) is the regular function 



Fe q (q 2 ) = 



L 2 2 + 2{^x/^r_ 

1 1 ] [l-(qx/u q ) 2 ] 2 



• 0J„ 



ArgTh^- - 1 



which is similar to the function F(q 2 ) found in the non-equilibrium case given by eq. (4.25). We can estimate 
the temp erat ure dependence of the photon density in the high temperature limit r»m: in this limit the 
integral (5.6) is dominated by momenta q ~ T and we can replace uj q and F eq (q 2 ) with their asymptotic 
expressions 



leading to the result 



q, F eq (q 2 ) 



shA z 

m T 







N v h{k, oo 



T^m^T^k 2e 2 T 3 

Tk 3 



In — + 0(1) 
m 



(5.7) 



B. Discussion 



Here we highlight a fundamental difference between our analysis of the photoproduction process and the 
typical analysis offered in the literature JIi"|-fl5[ . 

Our approach hinges upon computing the expectation value of the number operator of transverse photons 
in a state that has been evolved from an initial time to to the finite time t at which the number of photons 
is measured. By contrast, the usual approach computes the transition probability from a state prepared in 
the infinite past to a state in the infinite future. In such calculation there appears the familiar product of 
delta functions which are interpreted as the on-shell condition (energy momentum conservation) multiplied 
by the volume of space-time. Dividing by this volume one obtains the transition probability per unit volume 
and time which is interpreted as the production rate: this is basically the content of Fermi's Golden Rule. 
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In our approach we directly compute the expectation value (Nk)(t) = R^ + '(k,t) in a time evolved state 
and obtain the photon distribution at a time t by integrating this quantity, i.e., N p h(k, t) = f tg dt'R^ + \k, t'). 

This requires the knowledge of the dynamical photoproduction rate R^ + \k,t') for all times to < t' < t. 

The usual computation via Fermi's Golden Rule takes the long time limit and isolates the secular term 
that is linear in time by replacing R^ + \k, t') by its asymptotic limit 

R^(k) = lim R( + \k,t) (5.8) 



The condition (5.8) is tantamount to considering only on-shell processes, i.e, those that satisfy energy (and 
momentum) conservation. 

Keeping only on-shell processes, the large time limit of the photon number becomes 

N ph {k 1 t)=R ( +\k)-(t-t ) , £- t ^ oo. 

However our approach includes also off-shell processes that contribute to N p h(k,t) in a finite time interval. 
These processes do not contribute to N p h(k,t) asymptotically since they are subleading at very large time, 
i.e, 

t — >oo J J 

however they could be dominant at finite time. Actually, as we have seen in the previous section, the 
off-shell processes are of lower order in the electromagnetic coupling and strongly enhanced at soft momenta. 
Asymptotically we can write the the photon number in the form 

N ph (k, t) = N off - ahM (k, t) + R { +\k) -(i ~t ) 

where R^_ s h e n (k) is the usual rate calculated in equilibrium from on-shcll processes whose expansion in a 
begins at order a 2 (or aa s in the case of the quark-gluon plasma). In the case of broken symmetry studied 
in the previous sections 

N ff-aheu{k,t) cx a In \p(t - t )] 
Therefore off-shell processes dominate during a time scale t < t* with 

f *~J-lni 
/ia a 

This is an important point in the application of our novel approach to the physics of heavy ion collisions. 
In this case, the lifetime of the quark-gluon plasma is relatively short and the standard approach could miss 
important physics associated with transient off-shell effects. 

This analysis is essential in order to understand the possible phenomenological relevance of the transient 
effects. A quantitative assessment of it requires to compare the magnitude of the contributions to photon 
production from off-shell and on-shell processes at the finite time scale t of survival of the quark-gluon plasma. 
We intend to report the details of our studies on these issues within the context of photon production in the 
quark-gluon plasma in a forthcoming article. 

VI. THE MAGNETIC MASS OUT OF EQUILIBRIUM 

The magnetic mass in thermal equilibrium is defined as |n| 

ml qu ii, ma g = lim I™ (w) + KT/ ( 6 - 1) 

K — >U ^ — *U 
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where is the Fourier transform of the retarded transverse polarization kernel Yf^ ml {t — t') of the 

non-local part of the self-energy, and Sf'j* is the tadpole contribution in thermal equilibrium. When the 
evolution equation for the transverse mean-field is studied as an initial value problem, the relevant kernel to 
study is the Laplace transform of the retarded self-energy p8), i.e. 



ESS(*)= / dte-^^Zlit). (6.2) 



It is important to remark that the limits must be taken in eq.(S.l) in the precise order displayed above 
because the limits do not commute. 

It is a known result that in equilibrium the magnetic mass vanishes in an abelian gauge theory. The general 
argument relies on the structure of the Schwinger-Dyson equations, the Ward identities and translational 
invariance in space and time J30[ . More specifically to the scalar theory under consideration the vanishing 
of the magnetic mass to leading order in a (or alternatively to leading order in the hard thermal loop 
resummation) relies on the exact cancellation between the tadpole diagram and the zero frequency limit of 
the bubble diagram contributing to Yf^ utl {uj). A detailed analysis of this cancellation reveals the role of the 
Ward identity as highlighted by the general result in equilibrium. 

The equilibrium aspects of magnetic screening phenomena are fairly well established in abelian theories 
pO| , pl| , however to our knowledge the situation out of equilibrium has not received much attention. In this 
section we will study the dynamical aspects of the magnetic screening with an explicit computation at leading 
order in 1/N and first order in a. 

The initial stage in this program is to obtain an expression for the magnetic mass. 

This is achieved by considering the linearized evolution equation for the transverse photon condensate or 
mean field which is generated as a linear response to an externally prescribed transverse current Jt (x, t) . 
Such an equation has already been obtained in pS[ ] and we refer the reader to that article for details. In 
terms of spatial Fourier transforms, it is given by 



-^ + k 2 jA Tl (k,t) + J dt' B k ,ij(t,t') A Tj (k,t') = J Ti (k,t) , (6.3) 

where Y,k,ij{t,t') is the transverse retarded photon polarization out of thermal equilibrium. It contains two 
contributions, one local in time and determined by tadpole diagrams displayed in fig. and the other is 
non-local and retarded in time and given to lowest order in a by the bubble diagram displayed in fig. |^a. 
We have, 

= Zl ad (t) S(t - t') % + ^ b (t,t') Vijik), (6.4) 
where S^. a<i (t) is the tadpole diagram (fig. || b) 

Ej^t) = 2e 2 ($t$> = 2e 2 f A q 2 K(i)| 2 (6-5) 

Jo { 27r r 

and £^" b (£,£') the bubble diagram in real time (fig. || a), given by J|8| 

E»« 6 (f , f ) = -4e 2 J™ A q * dx {1 _ x 2 ) lm[G > (<) t , )G >^ (f } t/)] (6 6) 



with G < , G > the scalar Green's functions given by eqs.(3.58)-( |3.5S| ) and x being the cosine of the angle 
between q and k. 

The linear response to an external current is a different problem from that of photon production studied 
in the previous section, and although the polarization diagram shown in fig. |^ describes both processes, here 
we are interested in extracting a different information, which in equilibrium corresponds to the real part of 
the polarization in the limit of zero frequency. 
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Clearly, out of equilibri um, the very concept of mass is a delicate one, but we can make contact with the 
equilibrium definition |eq.( |6.l| )l by a derivative expansion in time. Writing 

= ^ ; r fc (t,f)= f x b k ub (t,t")dt" (6.7) 



and integrating by parts in eq.(6.3) we find 

(J^ + k 2 + 4 Qd (i) + r k (t,t)^ A Tl {k,t) - J dH r k (t,t') ^P(M') = J T i(k,t) 

Collecting the local terms in this equation of motion leads to the identification 



771 ^ 

mag 



limlim[Sf(t)+r t (M)] (6.8) 



We see from eqs.( p.l|) that t his d efinition reduces to the one at equilibrium in the case of time translational 
invariance. The definition (6.8) is the description of magnetic screening that is consistent with known 
equilibrium results in abelian theories. 

With the purpose of understanding the time and wavelength dependence of the several contributions, we 
now introduce a time and fc-dependent effective magnetic mass 

m 2 mag (k, r) = f dr' r') + E£*(r) = ml ub {k, r) + m 2 » (6.9) 

Jo 

The non-equilibrium definition of the magnetic mass, which coincides with the equilibrium definition in the 
case of time translational invariance is then 

m tnag = , lim , lim m L#. T ) 

and we remark again that the limits must be taken in this precise order. We now analyze generally the 
magnetic mass for both cases, broken and unbroken symmetry. 

We point out, however, that the effective magnetic mass ( |6.9| ) is i ntroduced to highlight the time scale 
of the different processes that contribute to the magnetic mass ( |6.8|) and its sole purpose is to provide a 
qualitative understanding of the different dynamical scales for the processes that contribute to magnetic 
screening. 

Using the asymptotic fo rm o f the mode functions (3.46)-(3.42), the definition of the asymptotic occupation 
numbers M q given by eq.( |3.60| ) and neglecting oscillatory terms that vanish in the asymptotic time regime 
due to dephasing, the tadpole contribution to the magnetic mass becomes 

2 rco 2 J 

m 2 tad (oo) = -^ 1M [1 + 2M] (610) 

2tt z ./n UJ„ 



which is reminiscent of the equilibrium tadpole contribution, but it contains the out of equilibrium distribu- 
tion functions Af q . 

The non-local contribution is given by 



Hub = lim lim T k (t,t) 



e 2 [°° f 1 r r 



(6.11) 



In the asymptotic time region, the mode functions are oscillatory and the product (pq(r) ip^ + ^{r) oscillates 
very fast for fcr> 1. Hence, any contribution that does not cancel the rapid time dependence of the phases 
will be averaged out. The memory integral from up to time r — > oo can be split into an integral from 
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t' = up to a time r' = tq > t^l within which the mode functions are exponentially growing but with slow 
oscillations and from r' = tq up to t' = r — ► oo. In this second region the mode functions have achieved 
their asymptotic forms (3.20) and ( [3.46 ). The contribution from the first domain cannot possibly cancel 
the fast oscillations from the mode functions at r. Therefore this first contribution will vanish by the rapid 
oscillation of the mode functions at very large r provided tr> 1. In the second r egion the integral can be 
performed using the asymptotic form of the mode functions and we find using eq.(3.6C) 



2 r>l. L fcr>l 
hub — 



e 2 .. ,. f°° q 4 dq f 1 2 
— - lim hm / — - — - — / dx (1 - ar) 

4tH k^or^ooj uj q J_ x 



{l - COS (u q + (t - t ) } 



UJ w+k\ - u i 



{l - COS [(« ]?f g| - WgJ (T - T )j } 



(6.12) 



This expression is remarkable, the terms with the occupation numbers are exactly of the same form as those 
obtained in an equilibrium description J28| and have a similar kinetic interpretation, the first describes the 
production minus the annihilation of two scalars, and the second term is out of equilibrium analogouos of 
Landau damping or bremsstr ahlun g (and its inverse) in the medium in terms of the asymptotic non-thermal 
occupation numbers A/" g (oo) ( 3.6C ). 

In the t — ► oo limit, the terms with cosines inside the integral ( 3.12| ) vanish. After taking the k — > limit, 
the integral over x is immediate and the Landau damping term leads to the derivative of the distribution 
function. Subtracting the vacuum contribution we find 



bub,ren 



e 



q 4 dq 



UJ q dAfq 

q dq 



(6.13) 



thus keeping k fixed and taking T->oowe recognize that the tadpole contribution and the one-loop bubble 
contribution have the same structure as the equilibrium calculation but in terms of the out of equilibrium 
distribution functions Af„ . 



Upon integration by parts in eq.( 6.13) a nd subtracting the vacu um contribution in ( 5.10 ) we find the exact 
cancellation between the bubble [eq.( |6.13 )] and tadpole [eq.(S.lO)] contributions, i.e. 



mag bub,ren 

just as is the case in the equilibrium calculation |2q 



tad^ren 



= 



(6.14) 



A. The effective magnetic mass m^ nag (k,T) for tr> 1 



Having established that the magnetic mass vanishes out-of-equilibrium, we can now study in detail the 
precise time evolution of the effective magnetic mass m 2 nag (fc,T) for late ti mes t ^> 1 and fixed but small fc, 
so that i;r» 1, by analyzing the differ ent co ntributions displayed in eq.(6.12). It is at this point that we 
justify keeping the oscillatory terms in (|6.12 ) so as to highlight the different time scales for the buildup of 
the different contributions. The two terms have very different oscillatory behavior, whereas the term with 
the sum of the frequencies maintains strong oscillations even if k 1, the second term proportional to 
the difference of the frequencies evolves slower in time for small k. This second term is recognized as the 
non-equilibrium analogous of Landau damping. In order to extract the long time behavior we proceed as 
follows: i) take k small and replace the difference in frequencies by a derivative with respect to momentum, 
ii) neglect the stron g osc illatory behavior arising from the term cos[2w (? t], to find that the effective magnetic 
mass defined by eq.(6.9) behaves as 



fer>l,T»l 



e 2 


[°°q 4 dq 1 


4tt 2 


Jo ^ q J 


e 2 


f°° q 2 dq 


7T 2 J 


'o w « 



' lu„ q dq 



1 — cos 



q kr 



(6.15) 
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where the first term inside the bracket is the contribution of the two-particle cut after neglecting the strong 
oscillatory component, the second term is the Landau damping term, the last term is the tadpole contribution, 
and we have subtracted as usual the vacuum contribution which is renormalized in the absence of the medium. 

Upon integrating over x and integrating by parts the Landau damping term, the time independent con- 
tributions cancel each other out as discussed above and eq. ( 5.15| ) yields 



fc«l,fcr»l 2e 2 



,, qdq ^{(krf 



and clearly in the long time limit, 



qkr\ qkr ( qkr 
sin I cos 



lim n»L( fc . T ) = 



qkr 



(6.16) 



in agreement with eq.(6.14). However this analysis clearly reveals that Landau damping or in-medium 
bremsstrahlung is the process with the slowest time scale in the long- wavel ength limit. 

The case of broken symmetry with u q = q is particularly clear. Eq.(6.16) then simplifies as 



m 2 maq {k,T) 



1 2 

/ct»i,t»i i e 



qdqAfq 



(k- 



rg [sin (fc r) — k r cos (k r)] 



(6.17) 



We see that m 2 nag (k,T) oscillates around zero for large kr with an amplitude that decreases as O f 
and period lix jk. 



B. The effective magnetic mass m^ nag (k,T) for fcr <C 1 



As we have noted above the asymptotic long time limit and the long-wavelengh limit do not commute, 
this happens out-of-equilibrium and also in equilibrium where the zero frequency and the zero momentum 
limit do not commute. 

However for finite time we can ask what is the behavior of the effective mass in the long- wavelength limit. 
This question is relevant for the evolution of the mean fiel d in the long-wavelength limit and for finite time. 
This corresponds to studying the effective magnetic mass (^1]) in the opposite limit kr <C 1 keeping t> 1. 

For finite time the effective magnetic mass is a slowly varying function of k thus for r 3> 1 but fcr < 1, 
we shall simply set k = to explore the region tel. 



In this case we find for the effective magnetic mass (6.9), 



m2 man( k i T ) 



fc<l,fer<l 



2 i-oo pT 2 rca 

q'dqj^ dr'Im[^(r)^ (rT ] + _^ f dq 



(6.18) 



For the computation of the bubble contribution rn 2 ub (k — 0, r) here [the first term in eq.( |6.18 )] some remarks 



are in order. Naively, since m 2 ub (k = 0, r) contains a product of four mode functions, each of order g~ 1//2 , 
one would expect the result being of order g~ 2 . However the following interference argument reveals that 
m bub(k — 0> r ) turns out to be of order 1/g, just as the tadpole contribution [the secon d ter m in eq.(6.18)]. 



Indeed, for < r' < tjvl the mode func tions fq(r') are given approximately by eq.(3.17). The exponen- 
tially gr owing term domi nates in eq.( 3.17 ) while the exponentially decreasing terms are of the order O(g) 
[sec cqs. (3.16) and (3.38)]. The dominant term has a time independent phase as can be seen read from 



Since the mode equations (3.11) have real coefficients, a solution with a constant phase during some time 
interval keeps such phase constant for all times. Therefore, the phase of the modes is time independent up 
to O(g) corrections and the phases of <fq( T ) an d l Pq( T ')* cancel up to O(g). Hence, Im^^T) ^(t')*] is a 
factor g smaller than g~ 2 . That is, it is of order g^ 1 and not of order g~ 2 . 
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Due to this cancellation of the dominant growing exponentials, an anal ytical evaluation of m 2 ^ ag (k,T 



requires detailed knowledge of O(g) corrections to the mode functions ( 3.17 ), which is not available analyt- 
ically. 



Instead, we evaluated numerically the integrals in eq.(6.18) using the high precision modes obtained in 
refs. ]25| , |37[ |. The results are displayed in figs. || and |7]. 

For late times, m^ aff (A: = 0, r) oscillates around the following constant values, 

e 2 

mL (k = 0, r 3> 1) = —0.000513 ... — for Broken symmetry, i]o = , 

9 

e 2 

m 2 (k = 0, t > 1) = —0.00126 ... — for Unbroken symmetry, tjq — 4 . (6.19) 

9 

2 

The coefficients of ^- are not very sensitive to the value of g for small coupling g <C 1, in the unbroken 
symmetry case the coefficient depends on the value of 770. 

Such small numbers arise from a delicate cancellation of the negative contribution from the bubble diagram 
and the positive contribution from the tadpole diagram. In the unbroken case, the larger is 770 the more 
negative is rr? mag (k = 0, r » 1). 

The negative sign of this effective squared mass indicates the unexpected presence of a weak instability 
in the time evolution of the mean field, which we conjecture to be linked to the strong photon production 
during this time scale. We expect to report on a detailed study of these issues in a forthcoming article. 

At this point it is important to remind the reader that had we studied a situation in which the global 
gauge symmetry was spontaneously broken either by the initial state or by the dynamics, there would have 
been a magnetic mass generated via the ordinary Higgs mechanism. Hence the vanishing of the properly 
defined magnetic mass is in agreement with the fact that the gauge symmetry is not spontaneously broken 
by the dynamics. 



VII. SCREENING AND DEBYE MASS GENERATION OUT OF EQUILIBRIUM 

The Debye mass or inverse of the electric screening length, determines the spatial extent over which electric 
charges are screened in the plasma. As in the case of the magnetic mass, the Debye mass can be obtained 
from a linear response problem. In this case the relevant linear response is that of the Lagrange multiplier 
Aq associated with the Coulomb interaction, or longitudinal photon to an external charge density ^Jo(x,t). 
The Debye mass can thus be recognized from the equation of motion for the expectation value of Aq(x, t) 
as a linear response to the external charge density. In terms of spatial Fourier transforms and calling this 
expectation value Ao(k 2 t) the equation of motion in linear response is obtained by following the method 
described in detail in [£8| . We obtain the following equation of motion for the expectation value in linear 
response 



k 2 A (k,t) + dt' ^(t,t') A (k,t') = J (k,t) , (7.1) 
Jo 

where the longitudinal retarded self-energy T,%(t, t') is given to lowest order in e 2 by the following expression 

n 

££(t, f ) = -4e 2 J J?L Im [d t ,G>(t, t') 9 t G> +g| (t, f ) - d t d v G> (t, f ) G> +] . { (t, t')] , (7.2) 

and j7o(fc,t) is the spatial Fourier transform of the external source that generates the linear response. 

We remark that Schwinger terms arising from the time derivatives of time ordered Green's functions had 
cancelled the tadpole contributi on 2 e 2 (<f>^<I>) and after this cancellation the remainder of the longitudinal 
photon polarization is given by (|7.2| ) • The reader is referred to J|^] for further details of this cancellation 
which is independent of whether the system is in or out of equilibrium. 
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Following the arguments presented previously in the case of the magnetic mass above, we define the Dcbye 
mass out of equilibrium as 



m 2 Deb = lim lim 

fc— >0 T— »0O 



f dr' ££(t,t') 
Jo 



(7.3) 



We emphasize again that the limits r — > oo and k — > must be taken in the order specified above since they 
do not commute. Taking the limits in the inverse order yields a vanishing result. 

Using the expressions for the Green's functions in terms of the mode functions as given by eqs. ( 3.5^ )-( [j~5"£ ) 
(in terms of the dimensionless mode functions Lp q (r)) we finally obtain 



m 2 Deh = lim lim e 2 



d 3 q 
(2^)3 



Im 



dr 



'Mr')> im (r'r\ (7.4) 



We compute now the Debye mass generally in both cases under consideration: broken and unbroken symme- 
try. The expression (7.4) displays a remarkable feature: unless the memory integral develops singularities as 
k — ► the Debye mass will vanish identically in this limit. Separating the time integral into a part between 
t' = and r 1 = tnl and a second part from t^l to r we recognize that no singularities can arise from the 
first part. The contribution to the Debye mass from the region of spinodal or parametric instabilities is reg- 
ular in the limit k — > and do not survive in the r — > oo limit. Therefore, we conclude that despite the fact 
that there are strong non-equilibrium processes during the stages of spinodal and parametric amplification, 
they are not directly associated with the generation of a Debye mass. However, as it will become evident 
below, the late time distribution of particles produced during these stages determines the Debye screening 
mass. In the second part of the integral the modes acquire their asymptotic form. Just as in the discussion 
of the magnetic mass, only few of the co ntri butions survive the rapid dephasing in the limit r — > oo. 

Replacing the mode functions in eq.(7.4) by their asymptotic behavior, using the relation ( 3.60 ) and 
neglecting the oscillatory contributions in the limit of r — > oo with k fixed, we obtain 



l Deb 



lim lim 

k— >0 t-^oo 



rf 3 
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(2tt) 3 w 9 



■JVc 



which is recognized as the longitudinal polarization evaluated at zero frequency 
contributions have an obvious kinetic interpretation which has been discussed in ref. 
Taking the zero momentum limit, we finally find 




b Deb 



djVg 

q UJq ~dq 



(7.5) 

Again the different 
(7.6) 



This expression reveals at once the important feature that the Debye mass is determined by the derivative 
of the distribution function of the charged fields with respect to momentum. Although this happens in other 
contexts, it is seldom highlighted in the literature. Integrating by parts at finite times the surface term at 
q = vanishes since the distribution Af q is regular at q = for finite times. We obtain the final form 



/??. 



Deb 



dq 

UJ a 



[2q 2 +M 2 (oo)]M q 



(7.7) 



We now study each case separately. 



A. Broken symmetry 

In the broken symmetry case with uj q — q ; A^(oo) = 0, the distribution function M q is of 0(1/ g) in the 
region < q < 1 and near the origin behaves as M q (r) *~ 1/q for all times including r — > oo. Fig. @ shows 
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g q M q {r) vs. q. Since lim^o q 2 N q = wc are justified in neglecting the surface term i n eq.( 7.6) and eq.(7.7) 



is valid even for r — ► oo. 
broken symmetry case 



Using equation ( |7.7| ), the relation ( 3.6C ) and the sum rule ( 3.25 ) we find for the 



m 2 D 



g tt z 



r[l + 0(g)} 



(7. 



B. Unbroken symmetry 

As discussed in detail in section 3.2, in the unbroken symmetry case the distribution function at times 
larger than t^l is dominated by the peak of the non-linear resonances. The distribution function continues 



to evolve at long times with two marked peaks in the region of non-linear resonances (3.4C) inside which the 
amplitudes A q (r) ; B q (r) and consequently the distribution function J\f q grows with a power law in time. 
The width of these non-linear resonant bands diminishes in time, the resonance near q w 770/V2 becomes 
subdominant and the resonance in the region < q < \J K\/t becomes the dominant one, the peak growing 
in amplitude and the width of the resonance diminishing as time evolves. Fig. ^ displays g q M q {r) vs. q for 
different times. 

Our extensive numerical calculations shows that for times r ^> tjvl the distribution function takes the 
scaling form 

A^(t) = G(g2 9 T) with G(0) = 0. (7.9) 
9 1 

The function G only depends on time and q through the combination x = q 2 T. We plot G(x) as a function 
of x in fig. [H| Notice that G{x) is of order 0(5°). 

At finite times the distribution function J\f q is finite at q — and neglecting the surface terms in the 



integration by parts leading to (7.7) is justified. 

Since the deriv ativ e of the distribution function is dominated by this peak, we find that even asymptotically 
the Debye mass (|7.6| ) continues to grow with time as \Jt. Fig. [n] displays the Debye mass as a function 
of time, it is clear from this figure that the tre nd i s that of monotonic increase as yfr as one obtains inserting 



the scaling form of the distribution function (7J)) into eq.(7.7) 



d q n« r / 2 d X s, n(U 

uj q 2gM 2 (oo) J a x 3 / 2 



Notice that the integral converges for r = 00 since G(oo) is finite [see fig. fto|| . 

The reason for this increase is that the distribution is dominated by the peak near q w which continues 
to evolve as a consequence of the non-linear resonance with the width ever decreasing in time and the peak 
continues to grow. 

In the infinite time limit the distribution will be peaked at zero momentum behaving as 

»r 1 \ G(oo) 

J\f q (r = 00) - 

9 q 

leadin g to a divergent Debye mass because A4 2 (oo) ^ and thus the behaviour at q = makes the integral 



in eq.( [7.7D to diverge. 

This divergence suggests that higher order contributions in the electromagnetic coupling must be taken 
into account and perhaps a resummation of higher order terms can lead to a finite Debye mass, but clearly 
this possibility requires a more detailed study which is beyond the scope of this article. 
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VIII. NON-EQUILIBRIUM TRANSVERSE CONDUCTIVITY 



Consider applying an external transverse electric field £T,ext{x,t) = —AT,ext(x,t) with Ar,ext(x, t) an 
external transverse vector potential. The induced transverse current is obtained in linear response by coupling 
the external vector potential to the current in the Lagrangian density £ — * C + Jt • Ar,ext- The transverse 
current induced by the external vector potential is obtained in linear response in terms of spatial Fourier 
transforms as 



(4(M)) = i dt' {jUk,t)4(-k,t')) ret A^ ext (k,f) 



(8.1) 



■2) 



where (J^(k,t)J^(—k,t')) ret is the retarded correlation function given by 

(4(fc,t)4(-fc,t'))r-e« = (4''(fc,t)J+ J '(-fc,t')> - {Jp\k,t)J T - J {-k,t')) 

and the symbols ± refer to the time branches along the CTP contour. Comparing with the expression 
for the photon polarization in the equation of motion (6.3) we recognize that the retarded current-current 
correlation function is given at lowest order (0(a)) by 



i{J l T {k,t).P T {-k,t')) ret = -£?*(*. V»{k) 

with T, b k ub (t,t') given by (fj). 

Introducing the non-equilibrium transverse conductivity as follows 

a i J{t,l/) = V ij {k)<T k {t,l/) = -V ij {k) f dt" ££" 6 (t,i") ; t > t' 



(8.3) 



(8.4) 



integrating by parts in eq.( |S.4| ) and neglecting surface terms we obtain the linear response relation 

(4(M)) = / dt'a l J(t,t') 4(M') (8-5) 



Although the definition of the conductivity (S.4) may not look familiar, it is straightforward to confirm 
that in the equilibrium case it leads to the usual relation between the conductivity and the polarization in 
equilibrium. 

In thermal equilibrium the polarization is a function of the tim e difference and the system has been in 
equilibrium from t = — oo. Thus, extending the lower limit in (8.4) to t" = — oo and writing 



^equil 
^k,bub 



(*"*") = 



+ OO 



i ^equil ( \ 



iuj{t—t") 



it is straightforward to find the spatial and temporal Fourier transform of the conductivity to be given by 



~ equil / \ 



a 



M = 



(8.6) 



which is the usual relationship between the bubble polarization and the equilibrium conductivity at lowest 
order in a. 

Since in the out of equilibrium case under consideration the initial state at time t = is the vacuum 
and the plasma is g ener ated during the stage of strong non-equilibrium evolution, the time integral in the 
conductivity kernel (8.4) has the initial time (t = 0) as the lower limit. 

The explicit expression for the conductivity at leading order in a follows from eq. ( |6.6|) and is given by 



Cfe(j) t') 



e 



/•OO p + l 

I q 4 dq / dx (1 — x 2 ) Im 



•7) 



It is difficult to compute explicitly the conductivity in the full range of the two time variables, however 
we can provide explicit formulae in the relevant regimes 1 < r' < r < t^l and when both time variables are 
in the asymptotic regime r > t' > tnl for fixed k. 
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A. 1 < t' < T < TNL 



In this time regime there are no fast oscillatory solutions and we can simply set k = to obtain an estimate 
for the long-wavelength limit of the conductivity. 

Broken symmetry: In this case the mode functions in this time regime are given by ( 3.17 )-( |3~i~8" ). 
Furthermore the modes that grow the most are those for q ss for which a non-relativistic approximation 
q <C 1 is reliable. The computation of the con ductivity proceeds in three steps, i) recognize the terms 
that contribute to the imaginary part in eq.(6.6), ii) carry out the integral in the variable t", iii) perform 



the integral in the q variable in the saddle point approximation using the non-relativistic approximation 
for the mode functions. In the region 1 < r' < r < tnl the dominant term of the conductivity in the 
long- wavelength limit is given by 



^«o(t,t') = — 



,2-r' 



T 3/2 r 5/2 



e 2T - 



15 7T 2 



[l + 0(k 2 r)] 



We remark that since the expression for the polarization (6.6) involves four mode functions one would naively 
conclude that the polarization and the conductivity would be ~ e 4r , however these terms are real and do 
not contribute to the imaginary part in eq.(6.6). Therefore at the end of the exponential growth of long- 
wavelength modes at t ~ r' ~ t^l the conductivity is of order 0{e 2 / g) and positive since the first term 
dominates over the second one in eq.(BJl) for r > r'. 



Unbroken symmetry: For this case the mode functions are given by eqs.( 3.32 )-( 3.36 ). The calculation 
of the conductivity in this time regime follows the same steps as in the broken symmetry case, with the 
difference in the third step being that the saddle point in the g-integral is at the maximum of the Floquet 
index q* given by eq. ( |3.39| ) . After some tedious but straightforward calculation we find 



OfcRjoO, t') 



t»l>i'>t'>1 



e 2 Vq 
768tt 3 / 2 



= 2S( J?0 )7 



(Vl + 1) 



1/4 



1 



q 9 n3/2 



i 



5„2 
4 '0 



i + o[ <zM,- 



(8.9) 



where £(770 ) = 4g(7 ?0 ) y/l + r% [l - 4g(ry ) + 0{q 2 { m ))} . 

Hence, besides some quantitative differences, this result is qualitatively similar to that in the broken 
symmetry case above with the same conclusion in the order of magnitude of the transverse conductivity in 
the long- wavelength limit at the time scale tnl. 



B. tnl < t' < t, k fixed 



In this regime we can use the asymptotic form of the mode functions generally in both cases, broken and 
unbroken symmetry. When we studied the magnetic mass in the previous section, we have noted that the 
long-wavelength limit does not commute with the long time limit. Thus, we will consider the long time limit 
but keeping k fixed. Furthermore, since the conductivity is a function of two time variables, we will consider 
r > t' and both arguments r and r' larger than tnl- In particular we will consider that phases involving 
the sum of the time arguments vary more rapidly and therefore dephase faster than those that depend on 
the difference of these time arguments, effectively deciding that r — r' is slower varying than r + r' but 
both arguments are in their asymptotic regime. The calculation proceeds along the same steps outlined in 
the case of the magnetic mass, the product ipq(r) V ? ^ + fe( T ) is strongly oscillatory in the asymptotic regime. 
Keeping only the oscillatory factors in the diffe rence of time arguments r — r', neglecting terms that oscillate 
much faster than these and using the relation (3.60), we find 



o-fc(r, t') 



^/^/> (1 -* 3) { (l+ ^ f+ * 



COS 



K + w i 9 -+fe|)( T - T ') 



K+ W |g+fc|) 
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CAfe|-^)- 



^\q + k\ -^?)( T - T ') 



This expression can be written in a more familiar form by introducing the representation 



(8.10) 



( W |9+fe| ± U l) 

which leads to the final relation 



= I du) 



j(r-r') 

~2uj 



Ifffel 



(T k (T,T') = J dui 



(8.11) 



.12) 



where £fc(u;) is the imaginary part of the Fourier transform of the asymptotic real-time, retarded polarization 
out of equilibrium, 



<7 4 dq 



W 1 u) \q+k\ 



dx{\- x 2 ) {(l + N q + A/j^) + w g + w lfl . +j?1 ) - 5(w - w, - w^) 

(AWI - -^) [^- w ? + w lfffcl)-^ + w 9- w l^l)]} ( 8 - 13 ) 



Hence, we find the remarkable result that in the asymptotic limit the conductivity and the polarization are 
related in a manner akin to that in equilibrium (|3.6|), i.e., 



CTfc(w) = 



.14) 



we stress however, that there is non-equilibrium information in this relationship because the distribution 
function of the produced particles is o ut of equil ibrium. 

The expression for the conductivity ( [B.1E )-( 8.14 ) simplifies considerably for small k and in the static limit 
u> = 0. Similarly to the Debye mass, to lowest order in a the conductivity is determinated by the derivative 
of the distribution function, 



5*(Q) 



fe^O 



47T 2 k 



00 d\r 

q dq 

k/2 



dq 



.15) 



I n th e broken symmetry case the k — > limit of the integral can be computed by part using the sum rule 
( 3.25 ) leading to 



CTfc(O) 



fc =0 e 2 m 2 R 



An 2 g k 



This result can be compared with the equilibrium conductivity at to — simply by replacing Af q by its 
thermal counterpart n q = (e q l T — l)^ 1 which yields 

~equil/ri\ k — >0 C T 

<v (°) = 



12 k 



This shows a qualitative comparison between the high temperature limit in thermal equilibrium and the 
small coupling limit out of equilibrium, which in dimensionful units reads 

T 
\m R \ 
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We emphasize that this comparison is only qualitative and should not be taken as a direct relation between 
the two cases since the non-equilibrium distributions are very far from thermal. 

Consider for instance the conductivity in the unbroken phase: our analysis in the previous sections shows 
that the distribution function J\f q is strongly enhanced at small q as a consequence of the non-linear resonance 
at q — 0, i.e, N q ~ 1/q 2 for q — > 0. This non-equilibrium effect leads to a logarithmic enhancement for long- 
wavelenths 



,,2 



. / n s fe->o e 2 (m R \ 
o*(0) ~ m R In i^—j , 

which has no analog in the equilibrium counterpart. 

We remark that the non-equilibrium conductivity [eqs.(B.13)-(8.14)] is finite for all k including k = at 
finite time. 

Only in the uj — limit, which corresponds to an integral up to infinite time, the conductivity has a 
divergent k — ► limit, such is the case for the equilibrium static conductivity a^ qull (0). Thus whereas the 
non-equilibrium conductivity is well behaved for long-wavelengths at any finite time, the long time limit 
will require a resummation of diagrams that must include the the width of the charged scalar particles. In 
the static limit the finite mean-free path of charged particles will provide a cutoff for long-wavelength (long 
distance) propagation and will lead to a finite long-wavelength static conductivity. 

Our analysis reveals i) the initial stages of build up of the conductivity through the formation of the non- 
equilibrium plasma, ii) an asymptotic description at long times in terms of the non-equilibrium distribution 
functions which is akin to the equilibrium description. 



IX. CONCLUSIONS, DISCUSSION AND FURTHER QUESTIONS 

In this article we have studied the formation of a plasma of charged particles from the strongly out of 
equilibrium processes of a) spinodal decomposition (or phase separation) and b) parametric amplification. 
The model, scalar QED with iV-charged scalar fields and a U(l) photon, docs not only provide an arena 
to study the questions of the formation of the plasma, electric and magnetic screening, photon production 
and conductivity out of equilibrium, but also is phenomenologically relevant both in cosmology and in heavy 
ion collisions as a description of the chiral phase transition out of equilibrium. In cosmology an important 
consequence of this study is the novel mechanism of generation of primordial magnetic fields at the time 
scale of the QCD phase transition, whereas in heavy ion collisions the mechanisms studied here can lead 
to strong photon production with non-equilibrium distributions that could be an important signature of 
non-equilibrium effects associated with the chiral phase transition. 

Spinodal decomposition describes the early stages during a quenched or supercooled second order phase 
transition and the dynamics is determined by the exponential growth of the fluctuations with wavevectors in 
the unstable band. Parametric amplification of quantum fluctuations occurs during the stage when the order 
parameter is oscillating around the minimum of the potential with large amplitude. In this situation there are 
resonances that amplify exponentially quantum fluctuations with wavevectors in the regions of parametric 
instability. In both cases, the explosive exponential growth of charged fluctuations lead to the formation 
of a non-equilibrium plasma and to photon production and the generation of electric and magnetic fields. 
These unstabilities are shut-off by the non-linear field interactions which are systematically and consistently 
treated in the large N limit |3(||37j . Thus we have combined the large N limit that allows the non-perturbative 
aspects of the formation of the plasma and a novel kinetic description of photon production to study many 
relevant electromagnetic properties of the non-equilibrium plasma. Our conclusions and further questions 
can be summarized as follows: 

• Photon production: We have obtained a novel kinetic equation to study photon production strongly 
out of equilibrium to lowest order in a and to leading order in the large N expansion. We find that 
at the end of the linear stage dominated by the exponential growth of instabilities in both cases, 
spinodal decomposition and parametric amplification, the photon distribution function is peaked at 
low momentum with a typical photon density of 0(a/X 2 ) with A the scalar self-coupling. In the case 
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of a quenched phase transition we find that electric and magnetic fields generated during the non- 
equilibrium stage are correlated on distances given by a dynamical correlation length £ « \A/I m fl| 
for times t < tjvl with |m^| the (renormalized) mass scale of the scalar fields. These mechanisms 
of photon production could be an important source of primordial magnetic fields in the early universe 
at a time scale of the chiral phase transition t w 10~ 5 seconds after the big bang and temperature 
scales O(100)Mev, however before coming to definite conclusions on the cosmological implications, two 
important issues must be studied further: i) whether a strongly supercooled (quenched) chiral phase 
transition can take place, given that the relaxation time scales for QCD are much shorter than the 
inverse of the expansion rate of the Universe during the transition, and ii) the kinetic equation used 
to study photon production neglected an initial population and therefore stimulated processes, these 
must be taken into account fully in the case of the cosmological phase transition since the Universe is 
radiation dominated at that stage. 

Perhaps, phenomenologically more relevant is the case of the chiral phase transition in Ultrarelativistic 
Heavy Ion Collisions, since it is quite possible that in this situation the phase transition occurs out 
of equilibrium In this case the photons produced through spinodal decomposition would have a 
non-equilibrium spectrum that could be a potential experimental signature . In the case of broken 
symmetry we have found that the presence of massless particles asymptotically, lead to collinear di- 
vergences in the bremsstrahlung contributions in the medium. These infrared divergences result in a 
logarithmic growth of the photon density asymptotically. This growth is also present in the equilibrium 
case and points out to a breakdown of the perturbative kinetic equation. A thorough understanding of 
the photon distribution function in this regime requires a consistent resummation using the dynamical 
renormalization group |46f| . 

• Magnetic screening mass: We have introduced a definition of the magnetic mass out of equilibrium 
which is the natural generalization of the equilibrium case. We find that the magnetic mass vanishes 
through a cancellation mechanism similar to that in the equilibrium case despite the fact that the 
asymptotic distribution functions are non-thermal. 

To highlight the different processes that contribute to the magnetic mass and their widely different 
time scales in the long-wavelength limit we have introduced an effective magnetic mass that coincides 
asymptotically with the proper definition of the magnetic mass. 

We find that the non-equilibrium generalization of Landau damping begins to compete with the con- 
tributions from two particle excitations and mean-field on time scales that in the long wavelength limit 
are far longer than those for these processes. This effective magnetic mass displays memory effects 
that correlates the spinodal or parametric particle production at early times with the dynamics at late 
times. We also find some unexpected weak long-wavelength instability in the time evolution of the 
mean transverse gauge field, which we conjecture to be related to the strong photoproduction during 
the early stages of spinodal or parametric instabilities. 

• Electric (Debye) screening mass: As in the case of the magnetic screening mass, we define the 
electric (Debye) screening mass out of equilibrium as the natural generalization of the equilibrium case. 

In the case of spinodal instabilities we find that the Debye mass is given by rn Deb = 8 |jti,r| 2 e 2 /A + 
O(A ). 

In the case of parametric amplification we find that the Debye mass diverges asymptotically as y/r 
with a coefficient of the order C(e 2 A _1 ). This result is a consequence of massive asymptotic states 
and the presence of non-linear resonances |37j] that result in a peak in the distribution function of 
the charged particles that moves towards zero momentum and whose width vanishes at long times. 



Since the Debye mass is determined by the derivative of the distribution function (7.6) in the case of 
massive particles a distribution which is singular at small momentum such as the one resulting from 
the non-linear resonances gives a divergent Debye mass. 

A divergent result to first order in a suggests that a resummation of electromagnetic corrections using 
for example the dynamical renormalization group Ew must be carried out. 
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• Transverse electric conductivity (Kubo): The transverse electric conductivity is an important 
transport coefficient which in the case of primordial magnetic fields limits the propagation and corre- 
lation of these fields and in the QGP it enters in the calculation of Ohmic energy losses in the plasma. 
We have obtained the non-equilibrium conductivity from Kubo's linear response out of equilibrium. 
The electric conductivity is a complicated function of two time variables and the wavelength. We 
solved in detail the build-up of conductivity during the early stages of formation of the non-equilibrium 
plasma as well as the asymptotically long time regime. The long-wavelength conductivity builds up 
exponentially because of the instabilities that lead to the formation of the plasma, and at the end of 
the stage dominated by linear instabilities it achieves a magnitude <7fc~o w |tor| e 2 /A. At asymptoti- 
cally long times we find that the conductivity has a similar structure to the equilibrium conductivity 
but with non-equilibrium distribution functions replacing the thermal ones. We find that the electric 
conductivity stays finite for all momenta including k — at finite times. 

Our study offers a novel view of the electromagnetic response of non-equilibrium plasmas in a model that 
allows to extract quantitative and qualitative information and that also bears phenomenological interest 
from the point of view of generating seeds of primordial magnetic fields at the chiral phase transition in the 
early universe and of describing non-equilibrium aspects of pion-photon dynamics in heavy ion collisions. 

APPENDIX A: KINETIC EQUATION FOR THE PHOTON DISTRIBUTION 

If the charged scalar fields were in equilibrium the rate of photon production would be determined by the 
imaginary part of the Fourier transform in frequency and momentum of the polarization depicted in fig. ^. 
The expression for photon production in equilibrium has been obtained in jll] Q for the case of the QGP. 

In the situation under consideration, strongly out of equilibrium, the polarization is not time translational 
invariant and the frequency representation is not available. The time evolution of the photon distribution 
function must be obtained from a kinetic equation. The validity of a kinetic description requires a wide 
separation of time scales between the time scale over which the photon distribution function changes and 
that of the phenomena that is strongly out of equilibrium. In the case under consideration the non-equilibrium 
evolution of the scalar fields result from spinodal and parametric ^instabilities and these occur on fast time 
scales of order |tojj| _1 ln(l/<7), we expect that the change in the photon distribution will occur on time scales 
that are longer by at least a factor 1/a. Hence under the assumption of weak electromagnetic coupling, the 
photon distribution function will evolve much slower than the non-equilibrium dynamics of the charged scalar 
field. Under these circumstances a kinetic description is valid. Furthermore since the charged scalar particles 
are far-off shell, a simple Boltzmann equation for the photon distribution function will miss the important 
off-shell effects associated with the non-equilibrium evolution of the scalar fields. This point becomes more 
important during the stage of spinodal instabilities when there is no meaning to on shell particles. 

In this Appendix we obtain the kinetic equation for the photon number from first principles (see also 

By using the simplest definition for the photon number or phase space distribution (for homogeneous 
systems) 

N vh (k,t) = (2tt) 3 ^-^- = £ («i( fc )«A(fc)> = —(A T (-k,t)-A T (k,t) + k 2 A T (-k,t)-A T (k,t)), (Al) 

A=l,2 

one extracts the time derivative of the distribution function as follows: 

N P h(k,t) = — <l T (t 1 ,-fe)A r (t 2 ,fc)+i T (t 2 ,-fc)A r (t 1 ,fc) 



+k 2 A T (t 1 , -k)A T (t 2 , k) + k 2 A T (t 2 ,k)A T (t 1 , k) > 



ti=t 2 
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2k dt 2 



d 2 

0^2 +k 2 ) < A T (h,-k) ■ A T (t 2 ,k) + A T (t2,-k) ■ A T (ti,k) > 



Therefore the photon density rate can be rewritten as 

»r / . . i d \(d 2 
N ^ = -2kdf 2 



■^ + k 2 )(g>{k-MM) + G<{k-t l M)) 



where Qij represent the exact photon propagator not to be confused with the free propagator Qij . 

In order to simplify this expression we need the Schwinger-Dyson equations for £/> and Q < . Including a 
mean field contribution 5il 2 (t) = 2e 2 ($^ (£)$(£)) in the Hamiltonian we have 



dt 2 

and 
dt 2 



+ k 2 )g? i (k;t 1) t 2 ) = -8n 2 (t)g? i (k;t 1) t 2 )+ /[n^(fc;ti,t)a< i (fc ; t ) t2)-n+ 7 :(fc;ti,^++(fc;t,t2)]^ 



+ k 2 )g< i (k;t 1 ,h)=-5&(t)g< i (k;t 1 ,h)+ [-n^-(fc;ti J *)^(fc;t > * 2 )+lC-(fc;ti,t)e-r(A;;t,*a)]*. 



Using the definitions 

U++(t, t'; k) = Il> m (t, t'- k)0(t - t') + nf m (t, t'- k)6(t' - t) 
n^T(t, t'; k) = Uf m (t, t'; k)6{t' - 1) + n< m (t, t'; k)9(t - t') 

n^(*, t' ; k) = n> m (t, t'; k), n+-(t, t' ; k) = n<jt, t' ; k) 

the photon production rate can be rewritten in the form 

N(k,t) = lsn 2 (t)^- [S>(M,ti) + S<(fc;Mi)] tl=t 



[n> m (fc; t, t') d t g<jk; t', t) - n< m (fc; t, t') d t g> t (k ; t', t)] dt'. 



(A2) 



This expression is exact but formal. To make progress, we consider the first order in a by replacing the 
full transverse photon propagator with its free field form but with non-equilibrium distribution functions 
and neglecting the electromagnetic contribution to the Green's functions of the charged scalar field. If there 
is an initial non-zero photon distribution N(k,to) the free Wightman functions read 

Qfi (k; t', t) = ± VS) {e-^ 4 '- 4 ) [1 + N(k, to)} + e^'^Nik, to)} , 



Qfi (k; t', t) = ± Vn (k) {e ifc (*'-*) [1 + N(k, t )] + e-^'^N(k, t )} 



By inserting this propagator in the first term (contribution from the mean field) in (A2) we see that the 
mean field does not contribute to the photon production to O(a). To first order in a the production rate 
obtains a rather simple form 



N(k,t) 



r« (k, t, t') (1 + N(k, t )) - r w (fc, t, t') N(k, t ) 



(!)/ 



dt' 



(A3) 
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with the time dependent rates 



il 1) (M,t') = -A- \v ij {k) nl (1 \k-t,t') e-*^-*') +v ij (k) n< (1) (fc;M') e ik ^ 



2k 



The transverse self-energies are given to lowest order by 

UlijW) = -4ie 2 | -0^ q Tl q Tj G>(t,t') Gj +E (*,0 
n£«(t,0 = -4ie 2 | ^ 9Tl G<(M') G< +£ (i,i') 



(A4) 



where in the large N limit the Green's functions G < > are given by ( 3.58D - fl3.59p leading to the final form of 
the self-energies to be used to lowest order in a and leading order in the large N limit, 



7*' (A) n>W(k;t,t') = ie "j^ i (1 - cos 2 (?) f q (t) m') f w+Sl (t) f; ?+ - kl (t') 



V^(k) n< 1 \k;t,t')=V i i{k) U^>(k;t',t) 



>(i), 



(A5) 



In the case in which the initial state is the photon vacuum, i.e. N(k,to) = the expression (A3) simplifies 
considerably. Upon integrating eq.(A3) in time up to time t, 



N(k,t) 



N(k,t') dt' 



one obtains two terms each one with a nested double time integral which can be written as a double integral 
up to the time t by inserting a theta function. Upon relabelling the time variables in one of the terms 
we obtain the expression ( |4.5| ), which is valid to lowest order in a and for vanishing initial population of 
photons. The photon production rate obtained in equilibrium in pd|-ff4T| also neglects the photon population 
in the initial state as well as the stimulated emission and the loss process. In these references the photons 
are assumed to escape from the medium without thermalizing and the photon production rate is valid to 
lowest order in a em and to all orders in the strong coupling constant. The kinetic equation (4.5) is precisely 



the non-equilibrium counterpart of the rate obtained in these references which is valid also to lowest order 
in the electromagnetic coupling and to leading order in the large N. 

This equation is clearly only approximate since it neglects the buildup of the population of photons. As the 
photon distribution increases in time due to photon production there will be stimulated photon production 
from the Bose enhancement factor resulting in enhanced photon production, but also processes in which 
a photon present in the plasma can decay into two charged scalars as well as change in the population 
with momentum k by bremsstrahlung or inverse bremsstrahlung in the medium. These latter processes 
will result in a depletion of population of photons and must be accounted for by a more complete kinetic 
description provided below whi ch m ust eventually be studied numerically. However, if the initial state is 
the photon vacuum, we expect ( |4.5| ) to be qualitatively correct for early and intermediate time scales, since 
Bose enhancement and the loss term from photon annihilation and scattering will depend on the number of 
photons present in the appropriate region of phase space. Hence first photons must be produced requiring an 
0(a) and then the stimulated and loss processes will take place requiring another power of a. Thus we expect 
that the early stages of photon production through spinod al de composition or parametric amplification will 
be described reliably with the simplified kinetic equation (4.5), the late stages will require the full kinetic 
equation described below which will involve a n ume rical study. 



The equation for the change of population (4.5) as well as the more general kinetic equation (A3) do 
not account for the change in the photon population, since N(k,to) is the population at the initial time. 
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Under the assumption of a wide separation of time scales, which relies on the weak coupling expansion in 
a a dynamical renormalization group analysis |39t| leads to a resummation of this kinetic equation by the 
replacement (to lowest order) N(k,to) — > N(k,t), thus leading to the final form of the (lowest order) kinetic 
equation 

N{k,t) = (k,t)[l + N{k, t)} - R a) (k, t) N(k, t) (A6) 
with the time dependent forward and backward rates given by 



R^\k,t) = f dt' T^\k,t,t') 

J t 



(A7) 



This resummation is akin to the Markovian approximation introduced in ]28||4"7| ] and is justified as a consistent 
expansion in the electromagnetic coupling. 

The resummation implied by this kinetic equation is based on a dynamical renormalization group analysis 
of kinetics [^9p^ j valid under the assumption of the separation of time scales, between the time scale of non- 
equilibrium processes of the scalar fields and that of the evolution of the photon distribution function, which 
is justified for small a. Physically the process that giv es rise to this kinetic equation is the following [^6|,|^],[59| : 



evolve the system from the initial time to with (A3) up to a time t\ at which the photon distribution has 
changed by a small amount of 0(a): 

N(k,U) = N(k,t ) + N(k,t')dt' (A8) 



with N{k,t r ) given by (AS) in ter ms o f the photon distribution at to- At this time t\ r eset the occupation 



number of the photon states to (A8) and evolve up to a further time ti using (A3) but now with the 



occupation number at time t\. The dynamical renormalization group establishes the equation that performs 



this operation infinitesimally and leads to ( |Aq ) [|39 46 1 . As explained in ref . [[39|j47| the coarse graining results 
from neglecting off-diagonal correlations of the type aa, in the time evolution of the density matrix. 

A similar resummation scheme is implied by the semiclassical Boltzmann equation, in which if the occu- 
pation numbers are treated in lowest order, the change is linear in time. Replacing the occupation numbers 
by the time dependent ones in the Boltzmann equation leads to a resummation and exponentiation of the 
time series jl7). However, as discussed in j|7l the Boltzmann equation assumes completed collisions that 
result in a coarse graining in time and neglects all of the transient effects and dynamics on short time scales. 

In particular, for the case considered in the previous section with vanishing photon occupation number in 
the initial state, the photon distribution function at a given time t is given by 

N(k,t) = f dh R^(k,h) e ~ P tl l(kM)dt2 (A9) 

J t 



with j(k,t) = R_ (k,t) — R < "l\k,t) being the total time dependent rate to 0(a). Clearly eq.(A9) provides 



+ 

a resummation of the perturbative series as is generally the case in any kinetic description wherein the rates 
are calculated perturbativcly. The early time behavior of the grow th of photon population is obtained from 
( |A9| ) by approximating j(k,t) ~ 0, leading to the expression ( |4.5| ). A more detailed estimate of the photon 
population including the reverse processes and depletion for a ge neral range of momenta k will undoubtedly 
require a numerical evaluation of the memory kernels in eq. ( |A9| ) , this is clearly a formidable task. 

If the time evolution were slow for all fields, we could write down a closed set of coupled kinetic equations 
for the distribution functions of photons and charged scalars. However, the strongly out of equilibrium 
evolution of the scalar fields and fast dynamics associated with the spinodal and parametric instabilities 
prevent such a kinetic description as there is no natural separation of time scales for the evolution of the 
scalar fields. The evolution of the scalar fields is therefore taken into account fully through the large N 



equations of motion and enters in the Green's functions that define the forward and backward rates ( A7) 
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